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SUMMARY 


This  document  serves  as  the  final  report  of  the  research  program,  ”  Acoustic 
Waves  in  Complicated  Geometries  and  Their  Interactions  with  Liquid-  Propellant 
Droplet  Combustion,”  sponsored  by  the  Air  Force  Office  of  Scientific  Research, 
Grant  No.  91-0171.  The  major  purpose  of  this  research  is  to  develop  a  com¬ 
prehensive  theoretical  analysis  within  which  multi-dimensional  acoustic  waves  in 
baffled  combustors  and  their  mutual  coupling  with  liquid-propellant  combustion 
can  be  properly  treated. 

During  this  report  period  from  March  14,  1991  to  May  14, 1993,  a  unified  linear 
and  nonlinear  acoustic  analysis  for  baffled  combustion  chambers  has  been  developed. 
The  linear  investigation  proceeds  by  solving  a  generalized  wave  equation  derived 
from  the  conservation  equation  for  a  two-phase  mixture,  utilizing  perturbation 
expansions  of  flow  variables.  Nonlinear  acoustic  study  of  baffled  chambers  is 
conducted  by  means  of  an  approximate  analysis  which  takes  into  account  second- 
and  third-order  nonlinear  gas  dynamics  and  combustion  response.  Results  suggest 
three  important  effects  of  baffles  on  transverse  modes  of  combustion  instability. 
First,  transverse  waves  can  be  longitudinalized  inside  baffle  compartments.  This 
may  decouple  combustion  response  from  oscillatory  motions  if  the  processes  near 
the  injector  face  are  sensitive  to  transverse  variations  in  pressure.  Second,  severe 
restriction  of  the  transverse  component  of  acoustic  velocity  fluctuations  is  observed 
in  baffle  compartments.  This  may  be  key  in  stabilizing  systems  with  velocity- 
sensitive  combustion.  Third,  the  frequency  of  oscillation  is  decreased  with  the 
addition  of  baffles,  in  agreement  with  experiments.  One  potentially  destabilizing 
effect  of  baffles  is  the  observed  pressure  concentration  at  the  injector  face.  Pressure 


1 


sensitive  processes  in  this  region  may  enhance  instability.  Limit  cycles  are  observed 
for  both  two-  and  three-dimensional  baffled  chambers.  Their  most  remarkable 
characteristic  is  high-frequency  modulation  of  the  limit  cycle  amplitude.  In 
addition,  the  conditions  for  the  existence  of  spinning  oscillations  are  treated  in 
depth. 

Some  of  the  results  obtained  from  this  research  have  led  to  the  following  two 
publications. 

1.  Yang,  V.,  Yoon,  M.W.,  and  Wicker,  J.M.,  ’’Nonlinear  Pressure  Oscillations  in 
Baffled  Combustion  Chambers,”  AIAA  Paper  93-0233,  presented  at  the  31st 
AIAA  Aerospace  Sciences  Meeting,  January  1993. 

2.  Yang,  V.,  Yoon,  M.W.,  and  Wicker,  J.M.,  ’’Linear  and  Nonlinear  Acoustic 
Waves  in  Liquid-Propellant  Rocket  Engines,”  presented  at  the  First  Interna¬ 
tional  Symposium  on  Liquid-Rocket  Combustion  Instability,  January  1993. 
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CHAPTER  1 

INTRODUCTION 


Combustion  instabilities  have  long  plagued  liquid-propellant  rocket  engine 
designers  and  developers.  The  essential  cause  of  this  problem  is  the  high  rate 
of  oscillatory  energy  release  within  a  volume  in  which  energy  losses  are  relatively 
small.  If  compensating  influences  acting  to  attenuate  the  oscillations  are  weak, 
then  unsteady  motions  in  a  combustor  may  give  rise  to  excessive  chamber  pressures 
and  heat  transfer  rates,  resulting  in  decreased  performance  at  best,  or  catastrophic 
failure  of  the  engine  at  worst.  Table  1  (assembled  from  Refs.  1-8)  presents  some 
practical  examples  in  which  pressure  oscillations  have  had  serious  impact  on  specific 
engine  development  programs.  As  this  table  shows,  high  frequency  instabilities 
have  remained  an  insidious  problem  throughout  the  entire  history  of  U.S.  rocket 
programs,  particularly  in  liquid  oxygen  (L0X)/RP-1  and  chemical  storable  systems. 

Because  instabilities  arise  from  sources  entirely  internal  to  the  system,  an 
external  observer  perceives  the  result  as  the  dynamical  behavior  of  a  self-excited 
system  capable  of  exciting  and  sustaining  oscillations  over  a  broad  range  of 
frequencies.  Typically,  the  oscillations  grow  out  of  a  reinforcement  between  the 
inherent  noise  of  the  combustion  process  and  the  numerous  characteristic  acoustic 
modes  of  the  feed-system/combustion  chamber  combination,  rather  than  from  any 
external  influence. 

Three  inherent  characteristics  of  the  combustion  chamber  environment  provide 
ideal  conditions  for  the  initiation  of  instability.  First,  combustion  chambers  are 


Table  1.  History  of  Combustion  Instabilities  in  Liquid-Propellant  Rocket  Motors 
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almost  entirely  closed.  This  condition  favors  excitation  of  unsteady  motions 
over  a  broad  range  of  frequencies.  Second,  the  internal  processes  tending  to 
attenuate  unsteady  motions  are  weak  relative  to  the  driving  energy,  thus  providing 
circumstances  that  allow  sustained  modes  of  instability.  Third,  and  perhaps  most 
important,  the  energy  required  to  drive  unsteady  motions  is  an  exceedingly  small 
fraction  of  the  heat  released  by  combustion.  In  typical  instances,  less  than  0.1 
percent  of  the  heat  release  is  sufficient  to  initiate  and  sustain  the  severe  pressure 
oscillation.9 

1.1  General  Features  of  Combustion  Instability  in  Liquid-Propellant 
Rocket  Engines 

Compared  with  solid  rocket  motors  and  other  liquid-fueled  propulsion  systems, 
the  combustion  processes  within  a  liquid  rocket  engine  have  several  characteristics 
that  produce  distinctive  differences  in  the  pressure  oscillations  observed.  First, 
except  for  highly  volatile  fuels,  such  as  hydrogen  and  methane,  most  propellants 
are  injected  into  the  chamber  as  a  spray  of  liquid  droplets  and  undergo  a  sequence 
of  atomization,  vaporization,  mixing,  and  combustion  processes.  Each  process  may 
serve  as  either  a  driving  or  damping  mechanism,  depending  on  its  sensitivity  to  the 
local  flow  oscillations.  Second,  the  combustion  chamber  is  directly  connected  to  the 
propellant  supply  system  through  the  injectors.  Interactions  between  the  unsteady 
motions  in  the  chamber  and  the  oscillatory  behavior  of  the  feedlines  may  affect 
the  stability  characteristics  of  the  engine.  Third,  the  occurrence  of  instabilities 
depends  intimately  on  the  operating  sequence  of  the  engine.  For  example,  the 
sudden  liquefaction  or  gasification  of  liquid  propellant  in  the  feedline  during  the 
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transient  start-up  and  shut-down  phases  may  trigger  instabilities;10  hence,  design 
criteria  must  make  provision  for  transient  as  well  as  steady-state  mechanisms  for 
initiating  combustion  instabilities.  Fourth,  the  propellant  also  frequently  serves  as 
the  coolant  for  a  regeneratively  cooled  system.  The  heat  transfer  and  other  engine 
performance  requirements  may  have  a  strong  impact  on  the  stability  characteristics 
through  their  influences  on  the  propellant  flow  conditions  in  the  injection  system.  A 
well  known  example  is  the  triggering  of  instabilities  in  the  J-2  engine  by  decreasing 
the  gaseous  hydrogen  injection  temperature.11 

Combustion  instabilities  are  a  particularly  complex  phenomenon  for  which 
several  modes  of  oscillations  have  been  commonly  observed.  They  are  often 
classified  according  to  their  spatial  structures  and  driving  mechanisms  as  low, 
intermediate,  and  high  frequency  instabilities.  Low  frequency  instabilities,  also 
known  as  ‘chugging’  or  POGO  instabilities,  received  much  attention  during  the 
1950s  and  1960s  as  a  serious  problem  with  the  Jupiter,  Thor,  Atlas,  and  Titan 
launch  vehicles.  The  underlying  mechanism  lies  in  the  coupling  between  the  fluid 
dynamics  in  the  combustion  chamber  and  the  propellant  supply  system.  Other 
effects  such  as  propellant  pump  cavitation  and  gas  entrapment  in  the  feedline 
may  also  contribute  to  driving  POGO  instabilities.  Characteristically,  this  mode 
corresponds  to  the  vibrations  of  a  Helmholtz  resonator,  with  a  frequency  range 
between  a  few  and  several  hundred  Hz.  Since  techniques  for  preventing  their 
occurrence  are  now  well  established,12,13  ‘chugging’  instabilities  are  no  longer  a 
serious  design  problem. 

Intermediate  frequency  instabilities,  also  known  as  ‘buzz’,  are  more  prevalent  in 
medium  size  engines  (25,000  to  250,000  N  thrust)  than  in  large  engines.  They  are 
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evidenced  by  a  growing  coherence  of  combustion  noise  at  a  particular  frequency. 
Two  sources  may  give  rise  to  their  occurrence.  The  first  is  associated  with  the 
interactions  between  the  unsteady  combustion  in  the  chamber  and  a  specific  portion 
of  the  propellant  feed  system.  The  second  arises  from  the  coupling  between  entropy 
waves  and  acoustic  motions  in  the  chamber.  It  is  well  established  that  for  bi¬ 
propellant  rocket  engines,  fluctuations  of  the  reactant  mixture  ratio  may  lead  to 
entropy  waves,  which  then  interact  with  the  nonuniform  mean  velocity  field  in  the 
nozzle  to  generate  acoustic  waves.  In  general,  the  ‘buzz’  instabilities  are  more  noisy 
than  damaging,  with  the  amplitude  of  pressure  oscillation  rarely  greater  than  5 
percent  of  the  mean  chamber  pressure,  although  in  some  cases  detrimental  high 
frequency  oscillations  can  be  triggered  by  this  mode.  Comprehensive  discussions  of 
‘buzz’  instabilities  and  entropy  waves  are  given  in  Refs.  14  and  15. 

High  frequency  instabilities,  often  called  ‘screeching’  or  ‘screaming,’  are  the 
most  common  and  vexatious,  and  are  characterized  by  reinforcing  the  interactions 
between  acoustic  oscillations  and  combustion  processes  inside  the  thrust  chamber. 
Depending  upon  the  response  of  combustion  processes  to  chamber  oscillations, 
energy  can  be  fed  into  acoustic  waves  such  that  their  amplitude  grows.  The  most 
destructive  result  of  these  large  amplitude  pressure  excursions  is  the  increased  heat 
transfer  to  the  chamber  walls  and  injector  face.  Burn-throughs  can  occur  in  a  few 
seconds  or  less,  causing  complete  failure  of  the  engine. 

This  type  of  instability  is  usually  characterized  by  well  defined  frequencies 
and  mode  shapes  that  correspond  closely  to  the  classical  acoustic  modes  of  the 
chamber.  Fundamental  modes  of  high  frequency  instabilities  are  thus  categorized 
as  longitudinal  and  transverse,  according  to  the  spatial  character  of  unsteady 


motions  within  the  combustion  chamber.  Longitudinal  modes  propagate  along  in 
the  axial  direction  of  the  chamber,  with  no  variation  in  the  transverse  plane  and 
are  usually  observed  in  chambers  with  large  length-to-diameter  ratios.  Transverse 
modes  exhibit  no  axial  variation  in  oscillatory  behavior,  but  propagate  radially  and 
tangentially  along  planes  perpendicular  to  the  chamber  axis.  For  contemporary 
rocket  engine  designs  in  which  the  chamber  aspect  ratio  and  nozzle  contraction  ratio 
are  relatively  small,  pure  transverse  modes  of  oscillations  dominate  because  of  the 
effective  damping  of  the  longitudinal  oscillations  by  the  nozzle  and  the  distribution 
of  the  combustion  along  the  chamber.  For  this  reason,  combined  modes  comprised 
of  the  superposition  of  longitudinal  and  transverse  modes  are  also  rarely  observed 
in  typical  liquid  rocket  engines.  Furthermore,  longitudinal  modes  are  usually  less 
destructive  than  transverse  modes,  so  theoretical  studies  such  as  this  one  should 
fully  address  transverse  modes  of  combustion  instability. 

Each  transverse  mode  may  exist  in  three  forms:  radial,  standing  tangential, 
and  spinning  tangential  waves.  The  differences  between  standing  and  spinning 
tangential  modes  can  be  clearly  seen  in  Fig.  1.1,  which  shows  the  particle  trajectories 
for  the  standing  and  the  spinning  modes  of  oscillations  in  a  simulated  F-l  engine 
combustion  chamber.16  Particles  move  back  and  forth  for  the  standing  mode  of 
oscillation,  but  have  an  epicyclic  trajectory  around  the  chamber  in  the  tangential 
direction  for  the  spinning  mode  of  oscillation.  The  standing  mode  particle  trajectory 
can  be  explained  by  the  fact  that  the  acoustic  velocity  varies  oppositely  in 
consecutive  cycles  with  respect  to  the  fixed  pressure  nodal  surface.  The  situation  is 
similar  for  the  spinning  modes  of  oscillation,  except  now  the  pressure  nodal  surfaces 
rotate  with  angular  frequency  corresponding  to  the  modal  frequency;  therefore, 


Standing  Tangential  Mode 


Spinning  Tangential  Mode 


Figure  1.1  Particle  Trajectories  for  Standing  and 
Spinning  Tangential  Modes 


particles  move  in  a  net  circular  fashion  as  time  elapses.  Spinning  wave  motions 
seem  to  be  usually  more  detrimental  because  of  their  effectiveness  in  agitating  gas 
molecules  in  transverse  directions,  thus  enhancing  heat  transfer  to  the  chamber 
walls.  However,  spinning  waves  may  also  be  accompanied  by  an  increase  in 
combustion  efficiency,  presumably  due  to  accelerated  mixing  processes.  This 
phenomenon  of  increased  efficiency  accompanied  by  decreased  stability  is  a  common 
trade-off  in  rocket  engine  designs,  and  the  ability  to  design  for  both  efficiency  and 
stability  simultaneously  represents  one  potential  payoff  for  instability  research. 

Severed  possible  mechanisms  have  been  proposed  that  may  be  responsible  for 


driving  the  spinning  mode  of  cornu  us  t  ion  instability.  Of  the  various  intermediate 
processes  occurring  during  the  combustion,  atomization,  vaporization,  droplet 
interaction,  mixing  of  the  vaporized  propellants  and  chemical  kinetics  are  the  most 
sensitive  processes  to  the  oscillations  of  velocity  and  pressure.  For  discussion 
purposes,  the  non-steady  effects  considered  here  can  be  divided  into  two  groups: 
effects  associated  with  atomization  and  vaporization,  and  effects  related  to  the 
mixing  process. 

The  atomization  and  jet  breakup  process,  and  its  relationship  to  spray 
formation  through  pressure,  temperature  and  velocity  perturbations,  can  affect 
the  energy  release  characteristics  of  the  gas  phase,  and  thus  the  stability  of  the 
combustor.  The  vaporization  process  that  is  directly  related  to  the  local  pressure, 
temperature,  and  velocity,  will  be  affected  by  oscillations  in  these  quantities. 
Furthermore,  there  can  be  mixture  ratio  gradients  in  the  vapor  because  of  the 
stratification  of  the  liquid  spray  in  the  liquid  propellant  injector.  If  the  transverse 
acoustic  field  is  imposed  on  such  a  spray,  the  vapor  will  be  displaced  relative  to  the 
droplet,  causing  mixture  ratio  oscillations  in  the  vicinity  of  each  vaporizing  droplet. 
Hence,  there  will  be  an  oscillation  in  the  burning  rate,  which  can  couple  with  the 
acoustic  field  to  produce  a  spinning  mode  of  combustion  instability. 

When  the  vaporization  rate  becomes  extremely  high,  it  is  possible  that  a  droplet 
is  heated  rapidly  through  its  critical  temperature.  With  droplet  shattering,  clouds 
of  the  very  fine  secondary  droplets  are  rapidly  gasified.  In  such  cases,  the  burning 
rate  could  be  controlled  by  the  rate  of  gas-phase  mixing.  Again,  this  oscillation 
of  the  burning  rate  is  coupled  with  the  transverse  acoustic  fields  to  generate  the 
spinning  mode  of  instability. 
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1.2  Mechanisms  of  Baffle  Operation 

To  provide  a  cure  for  combustion  instabilities,  or  to  at  least  limit  the  amplitudes 
to  acceptable  values,  passive  control  devices  (such  as  injector-face  baffles,  acoustic 
absorbers,  liners,  cavities,  and  quarter-wave  tubes)  can  be  used.  In  particular, 
baffles  can  provide  significant  stabilizing  effects  on  pressure  oscillations,  especially 
on  spinning  transverse  modes  of  instability  in  liquid  rocket  engines,  and  have  been 
widely  used  since  1954. 17  A  typical  configuration  consists  of  flat  plates  extending 
into  the  chamber  perpendicularly  from  the  injector  face,  arranged  in  a  radial  and/or 
circumferential  pattern,  as  shown  in  Fig.  1.2. 

Three  mechanisms  have  been  proposed  to  explain  successful  elimination  of 
instabilities  by  baffles.  These  are  (1)  the  modification  of  acoustic  properties,  such 
as  oscillation  frequency  and  waveform,  in  combustion  chambers;  (2)  the  restriction 
of  unsteady  motions  between  baffle  blades,  and  the  subsequent  protection  of  the 
sensitive  mechanisms  for  instabilities;  and  (3)  the  damping  of  oscillations  by  vortex 
shedding,  flow  separation,  and  viscous  dissipation.  Specific  examples  of  the  first 
mechanism  include  reduction  of  the  transverse  pressure  gradient  near  the  injector 
face  and  the  blocking  of  wave  paths  associated  with  certain  modes.  Spinning  modes, 
for  instance,  are  practically  unheard  of  in  baffled  combustion  chambers  due  to 
the  baffle  blade  boundary  conditions.  The  shielding  function  of  baffles  mentioned 
above  helps  to  minimize  coupling  between  parts  of  the  combustion  processes  from 
oscillatory  motions.  This  is  important  where  velocity  sensitivity  is  a  concern,  or 
where  triggering  of  an  instability  in  the  baffle  region  is  a  possibility.  The  third 
mechanism  of  baffle  operation-energy  damping-can  either  attenuate  the  oscillations 
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injector 


Figure  1.2  Schematic  of  a  Liquid-Propellant  Rocket 
Engine  with  Injector-Face  Baffles 

to  within  acceptable  levels,  or  extract  enough  energy  so  that  the  reinforcing  feedback 
between  combustion  and  chamber  acoustics  is  destroyed.  It  should  be  noted  that 
it  is  possible  that  more  than  one  mechanism  operates  at  the  same  time  in  a  baffled 
chamber,  but  some  may  dominate,  depending  on  the  specific  engine  parameters. 
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1.3  Research  Objective 

At  a  recent  JANNAF  workshop,18  it  was  suggested  that  acoustic  wave 
interaction  in  the  baffled  region  is  one  of  the  most  important  parts  in  the  study 
of  combustion  instability  for  liquid-propellant  rocket  engines.  Because  combustion 
instability  is  a  consequence  of  the  sensitivity  of  combustion  processes  to  ambient 
flow  oscillation,  any  realistic  investigation  into  the  engine  stability  behavior  must 
accommodate  a  thorough  treatment  of  the  wave  structure  within  the  baffled 
combustion  chamber.  However,  because  of  the  complex  nature  of  acoustic  flow 
in  the  baffled  region,  the  effects  of  baffles  on  wave  motions  have  never  been 
completely  understood  quantitatively.  At  present,  there  are  no  well-defined  criteria 
for  the  selection  of  baffle  configurations  that  will  lead  to  the  stable  operation 
of  an  engine.  As  a  result,  most  designs  in  use  today  are  based  on  experiences 
with  similar  combustor  configurations,  propellant  combinations,  and  operating 
conditions,  thereby  making  the  development  of  a  new  system  a  costly  trial-and-error 
process.19  Since  injector-face  baffles  provide  the  most  significant  stabilization  effects 
on  pressure  oscillations,  a  basic  understanding  of  the  oscillatory  flow  structure  in 
that  region  appears  to  be  a  prerequisite  in  treating  combustion  instabilities  in  baffled 
liquid  rocket  engines. 

The  primary  purpose  of  this  work  is  to  develop  a  theoretical  analysis  within 
which  multi-dimensional  acoustic  waves  in  a  baffled  combustion  chamber  can  be 
properly  treated.  Emphasis  is  placed  on  the  combustion  instability  behavior  of  a 
baffled  combustion  chamber  and  the  damping  mechanisms  of  the  baffle  as  described 
in  the  previous  section. 
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Two  other  focuses  in  this  study  are  the  examination  of  nonlinear  spinning 
wave  motions  in  the  unbaffled  chamber,  and  the  development  of  a  nonlinear 
baffle  analysis.  The  treatment  of  nonlinear  spinning  oscillations  is  crucial  due  to 
their  relevance  to  baffle  analysis  and  their  frequent  occurrence  in  practical  rocket 
systems.  Since  no  analytical  work  capable  of  treating  nonlinear  instability  in  baffled 
combustion  chambers  exists,  introduction  of  such  a  methodology  is  highly  desired 
and  necessary  if  questions  regarding  nonlinear  oscillations  in  physical  chambers  are 
to  be  addressed. 

1.4  Early  Works  on  Experimental  and  Theoretical  Investigation 
of  Baffles 

Several  experimental  investigations  have  furnished  results  which  lend  valuable 
insight  into  the  physical  mechanisms  of  baffles.  Experiments  using  non-reacting  or 
‘cold’  flow  allow  detailed  research  of  fluid  dynamic  phenomena.  Such  a  study  was 
conducted  by  Wieber,20  who  examined  the  acoustic  behavior  of  a  6-inch  diameter 
cylindrical  chamber  under  the  influence  of  various  baffle  patterns,  as  well  as  different 
injector  shapes,  chamber  lengths,  and  nozzle  contraction  angles.  Decay  coefficients 
of  pressure  oscillations  corresponding  to  the  first  transverse  mode  were  measured  for 
unbaffled  configurations  and  compared  with  internal  surface  area- to- volume  ratios 
of  the  chamber.  A  positive  correlation  was  discovered,  suggesting  that  viscous  wall 
losses  were  responsible  for  acoustic  wave  attenuation.  When  baffles  were  added, 
however,  acoustic  attenuation  increased,  but  no  correlation  with  the  internal  surface 
area  was  found.  A  more  complicated  damping  mechanism  than  simple  skin  friction 
alone  is  required  to  fully  explain  the  functions  of  baffles.  Wieber  further  observed 
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that  transverse  modes  were  not  only  damped  much  more  effectively  by  baffles  than 
were  longitudinal  modes,  but  that  the  reduction  in  frequency  from  the  unbaffled 
case  was  also  much  stronger,  presumably  due  to  relatively  longer  acoustic  particle 
paths.  Another  cold  flow  investigation  into  the  workings  of  baffles  was  initiated  by 
Torda  and  Patel;21  their  experimental  set-up  elucidated  vortical  motion  within  a 
baffle  cavity  as  a  result  of  transverse  flow  across  the  open  end  of  the  cavity.  The 
presence  of  circulating  flow  was  clearly  observed.  Viscous  damping  within  these 
vortices  are  thought  to  be  a  possible  mechanism  of  baffle  operation. 

Empirical  efforts  based  on  reacting  or  ‘hot’  flow  tests  can  yield  data  more 
readily  applicable  to  actual  rocket  systems.  Hannum  et  al.22  experimented  17 
different  baffle  designs  for  a  hydrogen-oxygen  engine  with  a  20,000  lbf  thrust.  In 
tests  with  at  least  a  2-inch  baffle,  stability  was  improved  from  the  unbaffled  chamber 
cases.  For  baffle  compartments  above  a  certain  size,*  longer  blade  lengths  were 
required  for  stability  as  the  cavity  dimension  was  increased.  Below  this  size,  all 
baffle  configurations  required  the  same  blade  length  for  stability,  independent  of 
compartment  measurement.  Interpretation  of  these  data  is  not  conclusive,  but 
one  possible  explanation  is  that  baffles  protect  the  sensitive  combustion  zone  from 
pressure  and  velocity  fluctuations.  Specifically,  a  certain  blade  length  corresponding 
to  the  extent  of  the  sensitive  combustion  zone  could  be  required  in  all  cases 
as  a  minimum  for  stability.  For  larger  baffle  cavities,  a  lor.smr  blade  would  be 
needed  to  impede  the  influence  of  the  main  chamber  oscillaHo  ,  >  .  acoustic  waves 

*  Compartment  size  was  defined  as  the  maximum  compartment  dimension 
measured  parallel  to  the  injector  face.  For  cases  of  noncongruent  cavities,  an  average 
weighted  according  to  the  number  of  each  cavity  style  was  used. 
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propagate  more  easily  through  larger  areas.  Although  no  relation  between  stability 
and  internal  surface  area  was  presented,  speculation  of  skin  friction  as  a  damping 
mechanism  is  prompted  by  some  results.  In  one  case,  a  2-inch  3-blade  radial  baffle 
became  unstable  when  it  was  burned  and  eroded  back  to  1.6  inches  in  length. 
An  egg-crate  configuration  was  stable  for  only  0.5  inch  blades,  perhaps  due  to 
the  greater  amount  of  surface  area  per  length  available  for  viscous  loss.  This 
conclusion  is  difficult  to  substantiate,  since  similar  behavior  is  not  consistently 
observed  throughout  the  data. 

Vincent  et  al.23  conducted  tests  similar  to  those  of  Ref.  22,  but  used  a  smaller 
engine  with  storable  propellants.  Trends  were  found  that  mirrored  the  findings  of 
Ref.  22  in  spite  of  the  differences  in  the  experimental  engine  itself.  A  sweeping 
summary  of  important  outcomes  from  several  experimental  baffle  programs  was 
completed  by  Hefner.24  Among  those  surveyed  were  the  Apollo  Service  Module 
Engine,  First  Stage  and  Second  Stage  Titan/Gemini,  Tran  stage,  and  the  Gemini 
Stability  Improvement  Program  (GEMSIP).  Specific  findings  are  too  numerous  to 
mention  here,  but  clearly  illustrate  the  labyrinthine  nature  of  baffle  data  reduction, 
as  well  as  the  dearth  of  fundamental  interpretive  understanding  regarding  baffles. 
Hence,  it  is  essential  to  formulate  a  theoretical  framework  suitable  for  systematic 
treatment  of  baffles  and  appropriate  deciphering  of  experimental  baffle  test  results. 

The  first  analytical  treatment  of  acoustic  motions  in  baffled  chambers  was 
undertaken  by  Reardon25  in  an  effort  to  explain  stability  data  from  the  GEMSIP 
program.  The  model  considered  three-dimensional  oscillations  in  the  main 
chamber  with  only  one-dimensional  longitudinal  oscillations  allowed  in  the  baffle 
compartments.  Combustion  was  assumed  to  be  concentrated  at  the  injector 
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face  and  modeled  with  a  time-lag  theory.26  Matching  of  acoustic  admittances  in 
both  regions  at  the  interface  provided  sufficient  conditions  for  determining  the 
linear  characteristics  of  the  entire  oscillatory  field.  This  approximation  yielded 
a  moderately  good  estimate  of  the  normal-mode  frequency  reduction  caused  by  the 
baffles.  However,  since  none  of  the  continuity  requirements  on  acoustic  pressure, 
velocity,  and  energy  flux  across  the  interface  were  fulfilled,  no  accurate  results 
of  acoustic  wave  structures  can  be  obtained.  Implementation  of  this  model  in 
predicting  engine  stability  is  subject  to  question. 

Oberg  et  al. 27,28  developed  a  more  elaborate  model  that  takes  into  account 
multi-dimensional  wave  motions  in  baffled  cavities.  The  homogeneous  Helmholtz 
equation  for  linear  acoustic  waves  in  the  entire  chamber  formed  the  backbone 
of  the  analysis.  The  effect  of  mean  flow  was  ignored,  with  the  combustion 
response  and  nozzle  damping  treated  as  boundary  conditions.  Separate  solutions 
were  constructed  using  Green’s  functions  for  baffle  compartments  and  the  main 
chamber.  The  acoustic  pressure  and  axial  velocity  in  each  region  were  then  matched 
at  the  interface  by  means  of  a  variational-iteration  technique  to  determine  the 
complex  wave  numbers  of  unsteady  motions.  Reasonable  results  were  obtained 
for  frequencies  in  two-dimensional  chambers,  but  no  solutions  for  cylindrical  and 
annular  chambers  were  reported  due  to  some  numerical  difficulties.  The  major 
limitation  of  this  theory  is  its  erroneous  prediction  of  engine  stability  behavior 
rendered  by  oversimplification  of  the  physical  processes.  Essentially,  this  model 
treats  an  inviscid  flowfield  with  concentrated  combustion  at  the  injector  face.  The 
resulting  solution  predicts  a  pressure  rise  in  the  forward  region  of  the  chamber, 
leading  to  a  destabilizing  influence  of  baffles  in  accordance  with  the  Rayleigh 


criterion.29  A  worsening  of  stability  with  increasing  baffle  length  is  obtained, 
rather  than  improvement  as  found  experimentally.  To  circumvent  this  problem, 
an  empirical  modification  that  allows  spatial  variation  of  the  combustion  response 
factor  along  the  injector  face  was  later  incorporated  in  the  analysis  to  match 
experimental  observations.30 

In  light  of  Oberg’s  findings,  it  is  evident  that  alteration  of  acoustic  wave 
structure  by  baffles  alone  can  not  explain  the  mechanisms  responsible  for  baffle 
damping.  The  stabilizing  influence  must  be  associated  with  other  means,  such  as 
viscous  dissipation,  vortex  shedding,  and  distributed  combustion  response,  which 
can  override  the  driving  effect  of  wave  alteration.  Baer  and  Mitchell31  established 
a  linear  acoustic  analysis,  with  emphasis  placed  on  the  fluid  dynamic  losses  along 
the  baffle  blades.  The  formulation  assumes  concentrated  combustion  at  the  injector 
face  within  the  confines  of  the  time-lag  theory,  and  solves  for  the  acoustic  velocity 
potential  using  an  eigenfunction-expansion  technique.  The  strong  flow  variations 
near  the  baffle  tips  are  treated  analytically  by  means  of  a  matched  asymptotic 
expansion  method.  A  turbulent  boundary  layer  is  then  ascertained  in  the  vicinity 
of  a  finite  thickness  baffle.  Although  results  considering  viscous  loss  within  this 
boundary  layer  yielded  realistic  acoustic  decay  rates,  the  conclusion  that  such 
dissipation  is  the  dominant  influence  of  baffles  may  be  unjustified.  For  instance, 
vortical  structures  induced  by  acoustic  motions  along  baffle  blades  and  convected 
downstream  into  the  main  chamber  may  play  an  important  role  in  dissipating  energy 
and  deserve  a  systematic  investigation.  Thorough  investigation  of  these  vortices 
using  information  gained  from  an  acoustic  study  of  baffled  chambers  is  done  by 
Wicker.32 
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1.5  Early  Works  on  Nonlinear  Combustion  Instability 

Spinning  transverse  modes  of  instability  have  very  often  been  seen  in  the  com¬ 
bustion  chambers  and  it  has  been  found  to  be  related  to  the  coupling  of  combustion 
processes  with  the  oscillations  in  flow  variables.  Several  analytical  and  experimental 
studies14,26’33’34  have  confirmed  that  the  pressure-sensitive  time  lag  concept  (n  —  r 
model)  can  be  applied  to  the  analysis  of  transverse  as  well  as  longitudinal  mode  of 
high  frequency  instability.  According  to  Reardon’s  experiment,35  it  is  not  sufficient 
to  explain  transverse  modes  of  instability  completely  by  the  pressure  sensitivity 
theory  only.  An  additional  combustion  rate  sensitivity  must  be  considered  for  the 
transverse  pressure  oscillations.  He  found  that  the  perturbations  of  radial  and 
tangential  velocity  components  have  a  strong  influence  on  the  combustion  process 
rates  in  a  manner  analogous  to  that  of  pressure  perturbations  for  the  n  —  r  model. 
Reardon  et  al.36  predicted  that  the  tangential  velocity  perturbation  had  a  strong 
destabilizing  effect  on  the  spinning  mode  but  had  no  effect  on  the  standing  mode  of 
instability.  Oberg  et  al.37  also  reached  this  conclusion-that  simple  pressure  coupling 
cannot  describe  the  spinning  wave  motions  in  the  baffled  combustion  chamber  and 
that  some  means  of  including  the  velocity  coupling  effect  are  necessary.  Several 
investigators38,39  addressed  the  importance  of  this  mechanism.  Thus,  a  conclusion 
can  be  drawn  that  velocity  coupling  as  well  as  pressure  coupling  plays  an  important 
role  in  the  spinning  tangential  mode  of  oscillations;  therefore,  the  effects  of  velocity 
perturbation  must  be  taken  into  account  separately  from  the  effect  of  pressure 
coupling  in  the  nonlinear  spinning  analysis. 

The  analysis  of  nonlinear  instability  behavior  in  baffled  combustion  chambers  is 
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another  part  of  this  work.  Although  studies  that  have  been  conducted  in  the  pursuit 
of  understanding  nonlinear  phenomena  have  not  specifically  addressed  baffled 
combustion  chambers,  valuable  building  blocks  have  been  laid  in  a  foundation 
used  in  the  current  study.  For  instance,  Sirignano  and  Crocco40  first  treated 
nonlinear  wave  motions  in  a  combustor  in  1964.  Using  a  perturbation  scheme  with 
characteristic  coordinates,  they  found  that  shock  waves  are  the  source  of  nonlinear 
losses  responsible  for  the  existence  of  the  limit  cycle.  Crocco  and  Mitchell41  treated 
the  pure  transverse  mode  in  an  annular  chamber  of  very  small  radial  extent  using 
a  nonlinear  perturbation  method  with  stretched  coordinates.  In  their  study,  the 
unsteady  combustion  processes  were  represented  by  the  sensitive  time-lag  (n  —  r) 
model. 

Zinn  and  co-workers42,43  took  a  different  approach  to  the  problem  of  nonlinear 
transverse  oscillations  in  a  liquid-propellant  rocket  engine.  Using  Galerkin’s 
method,  with  the  expression  of  nonlinear  solutions  as  an  expansion  of  classical 
acoustic  modes,  they  were  able  to  derive  a  system  of  second-order  ordinary 
differential  equations  governing  the  time- dependent  amplitude  of  each  acoustic 
mode.  This  system  was  then  solved  numerically  to  find  the  amplitudes  and  thus  to 
predict  the  existence  and  behavior  of  limit  cycles.  Owing  to  certain  assumptions, 
including  the  form  assumed  for  the  mean  flow  field,  the  extent  to  which  the 
results  may  be  valid  for  the  other  combustion  systems  could  not  be  addressed. 
Furthermore,  conclusions  were  based  on  numerical  calculations  of  special  problems, 
thereby  obscuring  possible  generality. 

Independent  of  the  works  cited  above,  the  approximate  analysis  for  nonlinear 
unsteady  motion  was  developed  by  Culick44  and  allows  expression  of  the  compli- 


cated  nonlinear  problem  in  terms  of  simpler  equations.  In  addition  to  using  the 
Galerkin’s  method  and  spatial  averaging,  the  method  of  time-averaging  developed 
by  Krylov  and  Bogoliubov45  was  incorporated  to  produce  a  system  of  first-order 
ordinary  differential  equations  for  the  time- varying  modal  amplitudes.  This  allowed 
for  a  more  convenient  analysis  than  the  second-order  equations  provided  by  spatial 
averaging  alone.  By  expanding  the  pressure  field  in  terms  of  normal  acoustic  modes 
of  the  chamber,  Culick  showed  that  the  nonlinear  behavior  could  be  equivalently 
represented  by  the  equations  for  a  system  of  nonlinear  oscillators.  Application  of 
this  work  to  motions  in  combustion  chambers  was  discussed  further  by  Culick,46 
with  special  emphasis  on  problems  related  to  solid-propellant  rocket  motors.  These 
general  ideas  have  recently  been  summarized  in  Yang  and  Culick,47  Kim,48  and 
Culick  and  Yang.9  This  methodology  is  an  invaluable  tool  for  studying  the  effect  of 
nonlinear  gas  dynamics  on  the  formation  of  limit  cycles,  and  is  an  important  part 
of  the  present  work. 

1.6  Thesis  Outline 

In  Chapter  2,  a  generalized  wave  equation  for  a  two-phase  mixture  is  derived 
utilizing  perturbation  expansions  of  flow  variables  to  study  the  oscillatory  flowfield 
in  a  baffled  combustor.  The  formulation  allows  for  acoustic  wave  motion, 
droplet  vaporization  and  combustion,  mean-flow/acoustics  coupling,  and  two-phase 
interaction. 

In  Chapter  3,  linear  acoustic  behavior  is  studied  for  two-dimensional  rectangu¬ 
lar  and  three-dimensional  cylindrical  baffled  combustion  chambers.  The  unsteady 
motions  are  expressed  as  a  synthesis  of  transverse  eigenfunctions  of  the  acoustic 
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wave  equation  in  the  baffle  compartments  and  the  main  chamber  separately.  The 
oscillatory  fields  in  these  two  regions  are  matched  at  the  interface  by  requiring 
continuity  of  acoustic  pressure  and  axial  velocity.  This  procedure  eventually  leads  to 
solutions  characterizing  the  unsteady  flow  structures  in  the  entire  chamber.  Strong 
effects  of  baffles  on  the  wave  forms  in  a  combustion  chamber  are  clearly  shown,  and 
the  effects  of  the  baffles  on  stability  are  explained  from  these  results. 

In  Chapter  4,  nonlinear  spinning  transverse  oscillations  in  a  cylindrical 
unbaffled  chamber  are  studied.  Velocity  effects  are  taken  into  consideration  through 
a  linear  combustion  response  function.  Using  this  approach,  linear  coupling  between 
modes  due  to  the  combustion  response  is  found  to  exist  and  to  play  an  important 
role  in  the  spinning  modes  of  instability.  The  conditions  for  the  existence  and  the 
stability  of  the  limit  cycle  are  also  examined  for  several  combinations  of  the  spinning 
modes. 

In  Chapter  5,  nonlinear  acoustic  analysis  is  performed  for  two-dimensional 
rectangular  and  three-dimensional  cylindrical  baffled  combustion  chambers.  The 
full  nonlinear  solution  for  the  baffled  combustors  is  formulated  in  a  manner 
analogous  to  Culick’s  approximate  method.  By  expressing  the  nonlinear  solution 
as  a  series  of  the  linear  modes  obtained  from  the  linear  acoustic  analysis,  time- 
varying  amplitudes  can  be  assigned  as  coefficients  in  the  series.  When  the  entire 
expansion  is  substituted  into  the  original  nonlinear  equation,  integration  over  the 
chamber  volume  yields  a  set  of  second  order  ordinary  differential  equations  in  terms 
of  the  time-varying  amplitudes.  The  solution  of  these  equations  gives  information 
regarding  nonlinear  phenomenon  such  as  limit  cycles. 
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CHAPTER  2 

THEORETICAL  FORMULATION 


Probably  the  most  important  fundamental  characteristic  of  combustion  insta¬ 
bilities  is  that  for  a  first  approximation  they  may  be  viewed  as  perturbations  of 
classical  acoustic  motions.  The  principal  nertur  -  lions  axe  due  to  the  combustion 
processes,  the  associated  mean  flow,  and  •  boundary  conditions  imposed  at  the 
inlet  and  exhaust  planes.  Treating  these  perturbations  within  the  framework  of 
classical  acoustics  has  been  successful  because  the  main  departures,  while  crucial  in 
defining  the  real  problems,  are  often  small  perturbations  in  some  sense.  Therefore,  it 
is  desirable  to  seek  the  theoretical  formulation  that,  when  all  perturbations  vanish, 
reduces  directly  to  a  representation  of  classical  wave  motion  in  an  enclosure.  The 
analysis,  therefore,  starts  with  the  establishment  of  a  generalized  wave  equation 
for  a  two-phase  mixture,  with  emphasis  placed  on  the  behavior  of  the  liquid  phase 
and  its  interaction  with  the  gas  flow,  matters  that  are  common  to  all  liquid-fueled 
systems. 

2.1  Conservation  Equations  for  a  Two-Phase  Flow 

Analysis  of  flows  in  liquid  rocket  engines  must  be  based  on  the  conservation 
equations  for  a  gas  containing  liquid  droplets.  A  proper  analysis  must  account  for 
the  differences  between  fuel  and  oxidizer  and  for  a  broad  range  of  sizes  of  liquid 
droplets.  However,  because  of  uncertainties  in  the  actual  flow  properties,  it  is 
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inappropriate  to  use  a  completely  general  formulation.  To  simplify  the  analysis,  a 
two-phase  mixture  in  the  combustion  chamber  is  assumed,  with  a  mass-averaged  gas 
comprising  all  gaseous  species,  including  inert  species,  reactants,  and  combustion 
products.  The  liquid  phase  is  treated  as  a  mass-averaged  fluid  comprising  fuel 
and/or  oxidizer  droplets.  The  complete  set  of  conservation  equations  thus  comprises 
those  for  the  gas  and  those  for  the  liquid;  then,  the  equations  are  combined  to  form 
a  set  that  governs  the  motions  of  a  single  medium.  This  idea  was  first  applied  to 
solid  propellant  rocket  systems  by  Culick  and  Yang,9  and  this  analysis  follows  their 
approach.  The  conservation  equations  for  this  two-phase  mixture  axe  given  below. 


Conservation  of  Mass  (gas) 

?£l+V-(p,  u,)  =  u,  (2.1) 

Conservation  of  Mass  (liquid) 

^  +  v  •  (f»u,)  =  (2.2) 

Conservation  of  Momentum 

d  ~ 

^  (PgUg  +  PI U,)  +  V  •  {pgUg u,  +  Pimm)  +  Vp  =  V  •  r  „  (2.3) 

Conservation  of  Energy 

d 

—  ( pgtgQ  +  pthl0)  +  V  •  (pgUgeg0  +  piUlhto)  +  V  •  (pug)  =  Q  +  V  •  q  (2.4) 


where  the  subscripts  g  and  /  refer  to  mass-averaged  quantities  for  gas  and  liquid 
phases,  respectively.  Evaporative  and  reactive  conversion  of  mass  from  liquid  to 
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the  gas  phase  occurs  at  the  rate  u>i  ( mass /sec- vol.).  The  viscous  stress  tensor 
is  symbolized  by  t  v  and  specific  total  internal  energy  of  the  gas  and  the  liquid 
stagnation  enthalpy  are  represented  by  ego  and  hi0,  respectively.  Furthermore,  Q  is 
the  rate  of  energy  released  by  homogeneous  reactions  in  the  gas  phase  (energy /sec- 
vol.),  and  q  represents  the  heat  flux  (conduction)  vector. 

With  simple  manipulations,  the  momentum  and  energy  equations  can  be 
written 

Pg-fc  +  P9U9  ■  Vu9  +  Vp  =  V  •  rv  +  Fj  -  (u9  -  u i)ui  (2.5) 


P9Cv-£  +  pgCvug  ■  VTS  +  pV  •  u9  =Q  -f  Qi  +  V  •  q  +  $ 

+  (hio  —  e9o)d>j  (^-®) 

+  u9  •  (Ug  -  u t)ui  +  (u /  -  Ug)  •  F / 

where  $  is  thp  dissipation  function.  The  force  of  interaction  between  the  gas  and 
liquid,  F /,  is  defined  as 

F,  =  -«[^+u,Vu,]  (2.7) 

and  the  heat  release  associated  with  chemical  reactions  and  heat  transfer  between 
the  phases,  Q/,  is 

Q'  =  -p'[^r+a>-Vh'  1  <2-8> 

The  enthalpy  of  the  liquid  hi  includes  the  heat  release  associated  with  the 
transformation  from  liquid  to  gas. 
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The  density  p  for  the  mixture  is  defined  as 


P  =  Pg  +  Pi  =  PgCL+Cm) 


(2.9) 


and  the  liquid  phase  to  the  gas  phase  density  ratio  is  given  as  Cm  =  pi/pg.  Mass- 
averaged  specific  heats  for  the  mixture  are  defined  in  the  usual  fashion.49,50 


Cv  = 


Cv  -t-  Cm  Ci 


c„  =  ^ 


Cp  4-  CmCi 


(2.10) 


l  +  cm  ~p  1  +Cm 

The  system  of  governing  equations  is  completed  by  introducing  the  perfect-gas  law 
for  the  mixture, 

p  =  pRTg  (2.11a) 


where  R  is  the  gas  constant  for  the  mixture 


R  —  CV~CV 


(2.116) 


Strictly,  Cm  must  be  treated  as  a  dependent  variable,  since  any  nonuniformity  in  the 
mean  flowfield,  or  unsteady  motions,  will  cause  the  liquid  droplets  to  slip  relative 
to  the  gas.  However,  if  the  droplets  are  dispersed  uniformly  in  the  gas,  Cm  is 
approximately  constant  throughout  the  chamber.  This  assumption  is  true  in  the 
solid  propellant  rocket,  but  not  as  valid  in  liquid-fueled  systems.  But  even  though. 
Cm  may  vary  significantly,  the  mass-averaged  thermodynamic  properties  are  not 
greatly  affected.  Based  on  this  assumption,  therefore,  the  momentum  and  energy 
equations  may  be  written  in  a  more  convenient  form  involving  the  mass-averaged 
properties  of  the  two-phase  mixture. 


•  Vu g  +  Vp  =  V  •  Tv  +  6Fi  +  6\\iu>i 


(2.12) 
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^  rdT 


pCv  [~Qj~  +  up  '  VTp]  +  PV  •  u g  —  Q  +  SQi  +  V  •  q  -(-  $  -f  ( hio  —  epo)u>/ 


(2.13) 


—  up  •  6uiuji  +  6ui  •  F/ 


where 


■  dSut 


(2.14) 


SFi  =  -pi  +  6ui  ■  Vup  +  Su i  ■  Vu /  +  u g  ■  V«5ui] 

rdShl  , 

SQi  =  -pi [~Q~  +  •  V6h,  +  6 U/  •  V(CtTg)  +  u9  •  (2.15) 

and  6\ii  =  u*  —  uff,  6h[  =  —  C{Tg.  The  energy  equation  (2.13)  is  combined  with 

Eqs.  (2.1)  and  (2.2)  to  give  the  equation  for  the  pressure: 


dp  „  „  R 

—  +  U g-Vp  +  7 pV  ■  Up  = 


Q  +  SQi  +  V  •  q  +  $ 

+  [(ft(0  -  ego)  -  Up  ■  (5u,]ci;,  +  6ut  •  F, 

-RTg[V-(pgug)+V-(piui)] 


(2.16) 


+  RTgUg  ■  Vp  +  pRTgV  •  Up 

The  chief  purpose  of  the  preceding  exercise  is  to  establish  the  forms  of  the 
equations  that  account  for  the  presence  of  liquid,  and  that  will  provide  a  good 
first  approximation  for  the  speed  of  sound  for  the  unperturbed  motions,  namely, 


a  =  yffRTg  = 


7 


1  +  Cm  Pg 


(2.17) 


This  formula  explicitly  shows  that  the  propagation  of  small  disturbances  is  governed 
by  the  elasticity  of  the  gas  (related  to  the  pressure),  and  by  the  inertia  of  the 
two-phase  mixture,  represented  by  the  factor  (1  +  Cm)pg.  Now  the  conservation 
equations  can  be  expressed  in  a  form  emphasizing  the  view  that  combustion 
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instabilities  are  unsteady  motions  best  regarded  as  perturbations  of  classical 
acoustics.  The  framework  for  the  analysis  is  based  on  the  sum  of  the  continuity 
equations  (2.1)  and  (2.2),  the  momentum  equation  (2.12),  and  the  energy  equation 
(2.16)  written  in  terms  of  the  pressure,  with  the  source  terms  in  a  general  form. 

Continuity  Equation 

|+us.V„=W 

Momentum  Equation 

d\i„ 

P-Qf  +  Pug  ‘  -  -Vp  +  T 

Energy  Equation 

— ^  +  7pV  •  u g  =  -uff  •  Vp  +  V 
For  the  circumstances  treated  above, 

W  =  -pV  •  us  -  V  •  (piSui) 

T  =  V  •  +  <5F<  +  8x1^1  (2.22) 


(2.21) 


(2.18) 


(2.19) 


(2.20) 


^  R 
"  cv 


Q  +  8Q,  +  V  q  +  $+  [(ft/o  -  eff0)  -  u5  •  6u,]u>, 


+  <5uj  •  Fj  —  CvTgV  •  (pj6uj) 


(2.23) 
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2.2  Nonlinear  Wave  Equation 


To  derive  a  nonlinear  wave  equation  governing  the  unsteady  motions  in  the 
chamber,  the  dependent  variables  are  first  decomposed  into  mean  and  fluctuating 
parts, 

p(r,t)  =  p(r)-f  p'(r,<) 


ug(r,t)  =  uy(r)  +  u'g(r,t) 
p(r,t)  =  p(r)  +  p'(r,t) 


(2.24) 


It  is  assumed  that  the  mean  quantities  are  time-invariant  and  the  fluctuating  parts 
are  small  perturbations.  Although  variation  of  the  mean  quantities  is  typically  small 
in  actual  rocket  engines,  spatial  dependence  is  allowed  for  in  the  general  formulation. 

To  the  second  order  in  fluctuations,  Eqs.  (2.11a),  (2.18),  (2.19)  and  (2.20)  can 
be  rewritten  in  their  perturbed  forms,  giving  a  complete  set  of  equations  for  the  six 
unknowns:  p',  p',  T'g ,  and  the  three  velocity  components. 

^  +  us  •  vp'  +  u;  •  Vp  +  u;  •  vp'  =  w'  (2.25) 


(2.27) 
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P  =  R(p'Tg  +  pTg)  +  Rp'T'g  (2.28) 

Since  acoustic  waves  manifest  themselves  by  the  presence  of  pressure  oscillations, 
and  pressure  signals  are  easily  measured  and  processed,  the  wave  equation  can  be 
most  conveniently  written  in  terms  of  pressure  fluctuation  p'.  Consequently,  taking 
the  time  derivative  of  Eq.  (2.27)  and  substituting  Eq.  (2.26)  for  du'g/dt  yields 

v2p'-?l£  =  fc  =  '‘-  +  '‘’  +  '‘’  <229> 

where 


hi  =  -V  -  p(us  •  Vu'  +  u'g  •  Vu})  +  p'ug  •  Vus 


1  dp'  du'  dp' 

+¥  ^V  u’  +  -af-v?+a’-vW 


ho  =  -V 


+  : 


Q^l  >1 

•  Vuj)  +  P’*9  ■  Vu*  +  p'-fif 


|«-Vp')  +  7|(P'V'<) 


dt' 


(2.30a) 


(2.306) 


h3  = 


V  ■?'  - 


1  dV' 
a 2  dt 


(2.30c) 


Here,  the  subscripts  1,  2  and  3  represent  the  linear  effects,  second  order  nonlinear  gas 
dynamics,  and  the  effect  of  source  terms  such  as  combustion  and  viscous  dissipation, 
respectively. 

Boundary  conditions  set  on  the  gradient  of  p1  are  found  by  taking  the  scalar 
product  of  the  outward  normal  vector  with  the  perturbed  momentum  equation,  and 
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then  applying  appropriate  acoustic  admittance  functions  along  the  surface  of  the 
field. 


n  •  Vp'  =  -n  •  f  =  -n  •  (fj  +  f2  +  fa) 


(2.31) 


where 


r  <9u'  ■ 

fl  =  f  dt  +  p[n9  ■  Vus  +  us  • 

(2.32a) 

,du'a  ,  , 

f2=lpV +^ui-Vu;} 

(2.32  b) 

&* 

H 

i 

Si 

(2.32  c) 

If  all  perturbations  are  absent,  functions  h  and  f  vanish,  and  the  wave  equation 
for  the  pressure  in  classical  acoustics  with  the  boundary  condition  for  a  rigid  wall, 
n  •  Vp'  =  0,  is  recovered. 
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CHAPTER  3 

LINEAR  ACOUSTIC  WAVES  IN  BAFFLED  CHAMBERS 


The  effects  of  baffles  on  acoustic  motions  in  combustion  chambers  have  not 
been  understood  quantitatively  due  to  the  complicated  chamber  geometry  and 
associated  difficulty  of  analyzing  the  flowfield  in  the  baffle  region.  Since  the 
elimination  of  combustion  instabilities  is  the  most  important  but  least  understood 
aspect  of  baffle  design,  this  chapter  deals  with  the  linear  stability  behavior  of 
baffled  combustion  chambers  as  well  as  the  actual  physical  processes  responsible 
for  suppressing  unsteady  motions.  The  wave  equation  derived  in  Chapter  2  is 
applied  to  the  baffled  combustion  chamber  domain  to  investigate  the  acoustic  wave 
characteristics  for  both  two-  and  three-dimensional  chambers. 

3.1  Linear  Acoustic  Oscillations 

For  a  three-dimensional  cylindrical  baffled  combustion  chamber  as  shown  in 
Fig.  3.1,  a  direct  treatment  of  the  wave  equation,  subject  to  the  appropriate 
boundary  conditions,  appears  to  be  formidable  because  of  the  geometric  complexity 
imposed  by  the  baffles.  To  facilitate  the  formulation,  the  form  of  the  acoustic  field 
solution  is  constructed  in  two  parts:  the  baffle  compartments  and  the  main  chamber. 
Since  the  spatial  variations  of  acoustic  motions  are  quite  different  in  the  baffle 
compartments  and  the  main  chamber,  the  oscillatory  fields  in  these  two  regions  are 
best  treated  separately.  The  eigenfunction  expansion  technique  is  adopted  to  solve 
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Figure  3.1  Schematic  of  a  Three-Dimensional  Baffled 
Combustion  Chamber 

the  wave  equation,  since  the  standard  separation  of  variables  method  cannot  be 
used  due  to  the  boundary  shape  of  the  domain. 

The  fluctuating  pressure  in  each  region  is  synthesized  as  a  Fourier-type 
expansion  in  terms  of  eigenfunctions  for  the  cross  section,  but  allows  for  temporal 
and  axial  variations  through  the  series  coefficients. 

OO 

P'(M)  =  P  53CnOM)V>„(r,0)  (3.1) 

n=0 

where  the  eigenfunction  xj)n  satisfies  the  Helmholtz  equation  in  the  transverse  plane: 


=  0 


(3.2) 
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subject  to  the  boundary  condition 


o 

II 

c 

l> 

a 

(3.3) 

The  transverse  Laplacian  operators  V2  in  the  Cartesian  and  cylindrical  coordinates 

are  respectively  defined 

as 

_2  d 2  d 2 

Vr  “  dy 2  +  dz 2 

(3.4a) 

and 

^=ld(rd\l  S2 

T  r  dr  \  dr  J  r2  d6 2 

(3.46) 

Then,  the  complete  Laplacian  operator  V2  becomes 

V2  -  (3  5> 


The  subscript  n  stands  for  the  doublet  of  indices  (m,  s)  defined  below.  Due  to 
the  azimuthal  degeneracy  in  the  cylindrical  chamber,  has  two  possible  forms: 

V’n  ~  ^  ^  (3-6) 

t  sin  mo 

The  transverse  eigenvalues,  represented  by  the  wave  number  /cma,  axe  determined 
by  the  roots  of  the  derivative  of  the  mth  Bessel  function  at  the  chamber  wall.  That 
is, 

dr 

Now  multiply  the  wave  equation  (2.29)  by  il>n  and  Eq.  (3.2)  by  p',  subtract  the 
results,  and  integrate  over  the  cross  section  to  find 

ds  =  j J 'pnhds  +  /  ^  • ndl  (3-8) 
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where  the  line  integration  on  the  right  hand  side  is  conducted  along  the  baffle  blades 
and/or  the  circumference  of  the  chamber  wall  at  each  cross  section,  and  the  function 
fT  is  defined  in  the  transverse  plane  as 

fT  •  n  =  VTp'  ■  n  (3.9) 

Substitute  Eq.  (3.1)  into  Eq.  (3.8)  and  rearrange  the  result  to  get  the  following 
differential  equation  for  £„(x,  t). 

-  «2| Jr  +  =  -Jei  [J I  ** hdS  +  /  • ndl  (3-10) 

where  a is  the  frequency  of  the  normal  transverse  mode  ipn ,  defined  as 

wn  =  (3.11) 

The  Euclidean  norm  of  the  mode  function,  En  is  defined  as 

El  =  jji!>ldS  (3.12) 

For  linear  acoustic  problems,  all  the  property  functions  and  variables  are 
assumed  to  vary  in  a  time-harmonic  fashion, 

P'(r,t)  =  p(r)e,nt 
u;(r,t)  =  us(r)ei«‘ 

Cn(x,<)  =  Cn(x)e‘m  (3.13) 

h(r,t)  =  h(  r)e,n< 
fT(r,8,t)  =  f r(r,0)e*nt 
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where  Q,  is  the  complex  frequency  of  oscillation  and  is  defined  as 


Q,  =  u>  —  ia 


(3.14) 


The  primary  purpose  of  linear  stability  analysis,  then,  is  to  determine  the  radian 
frequency  u>  and  the  growth  constant  a. 

The  ordinary  differential  equation  for  the  axial  distribution  function  is 
obtained  from  the  substitution  of  Eq.  (3.13)  into  Eq.  (3.10). 


d2  Cn 

dx 2 


+  ^nCn  =  ~ 


PEI 


J J  IpnhdS  +  £ 


V'nfr 


ndl 


(3.15) 


where  the  modified  wave  number  6n  is  defined  as 


“  a2 


(3.16) 


Combustion  is  often  assumed  to  be  concentrated  at  the  injector  face,  an 
assumption  that  is  based  on  the  experimental  observation  that  the  majority  of 
the  combustion  processes  are  completed  near  the  injector  of  the  chamber.51  Given 
this  condition,  the  mean  flow  Mach  number  is  considered  as  constant.  That  is, 


us  =  Maez  (3.17) 

where  M  is  the  constant  Mach  number  and  ex  is  the  unit  vector  in  the  axial 
direction.  However,  this  assumption  consequently  limits  the  range  of  mean  flow 
Mach  number.  For  example,  if  the  Mach  number  of  the  average  flow  is  larger  than 
roughly  0.4,  the  approximation  of  classical  acoustics  deteriorates  since  the  Doppler 
effect  and  refraction  may  cause  substantial  distortions  of  the  acoustic  field.  But  the 


mean  flow  Mach  number  is  generally  small  enough  over  most  physical  combustion 
chambers.  With  this  assumption  the  linear  mean  flow  effect  hi  in  the  wave  equation 
can  be  expressed  as 


hi 


-pv 


u 


du'g 
1  dx 


1 

+ 

a 1 


9 


(3.18) 


Within  first  order  accuracy,  the  acoustic  velocity  in  Eq.  (3.18)  can  be  expressed 
in  terms  of  the  acoustic  pressure  p'  using  the  classical  acoustic  momentum  equation. 


du’g 

dt 


(3.19) 


Substitution  of  Eq.  (3.19)  into  Eq.  (3.18)  yields  the  final  expression  for  the  linear 
effects  of  the  wave  equation. 


hun  ~  - 


iug  d  /X72-\  .  r-,- 

-£tx(Vp)+  ¥u’Vp 


(3.20) 


Substitution  of  expression  (3.20)  into  Eq.  (3.15)  leads  to  a  second-order  ordinary 
differential  equation  for  the  axial  distribution  function  Cn* 

2MSlid(„  ,  „2;  „ 

'^“'T_'57+S"Cn  =  0 


The  major  task  at  this  point  is  to  determine  the  complex  frequency  character¬ 
izing  linear  pressure  oscillations.  Two  cases  are  treated  separately:  two-  and  three- 
dimensional  baffled  chambers.  For  both  cases,  the  acoustic  wave  characteristics  are 
investigated  systematically  to  give  a  basic  understanding  of  the  wave  structures  in 
the  chambers. 
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3.2  Two-Dimensional  Baffled  Chamber 


Although  the  configurations  of  most  practical  baffled  combustion  chambers 
are  three-dimensional,  it  is  beneficial  to  first  achieve  the  basic  understanding  of  the 
oscillatory  fields  in  a  two-dimensional  chamber.  This  model  simplifies  the  geometric 
complexity  of  the  problem,  and  clearly  presents  some  important  mechanisms  of 
baffles  for  modifying  the  characteristics  of  unsteady  motions.  Fig.  3.2  shows 
the  schematic  of  a  two-dimensional  baffled  combustion  chamber.  The  acoustic 
fields  in  the  baffle  compartments  and  the  main  chamber  are  linked  together  by 
requiring  continuity  of  acoustic  pressure  and  axial  velocity  at  the  interface.  This 
procedure  eventually  leads  to  a  transcendental  equation  for  the  complex  wave 
number.  Results  provide  explicit  information  about  the  acoustic  wave  structure 
and  complex  frequency. 

3.2.1  Acoustic  Field  in  Baffle  Compartments 

Because  Eq.  (3.21)  is  a  second-order  ordinary  differential  equation  with 
constant  coefficients,  the  solution  in  a  two-dimensional  model  takes  the  form 


where 


Cmb(x)  =  Ambtqmbl  x  +  (3.22) 


9m  k,l  - H 

a 


Mfi 


a 


T  2 


+  oi 


mb 


MQ  _ 

9m,, 2  - 


MQ 


a 


+  01 


mb 


(3.23a) 

(3.236) 
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Figure  3.2  Schematic  of  a  Two-Dimensional  Baffled 
Combustion  Chamber 


The  subscript  b  stands  for  baffle  compartments.  The  acoustic  boundary  condition 
at  the  injector  face  is  characterized  by  an  acoustic  admittance  defined  as 


A/==i!Z L 

P'hP 


(3.24) 


x=0 


This  leads  to  a  relationship  between  the  coefficients  Am(,  and  i?m&. 

r  _  Bmb  _  9m»,»  +  + 

Lsmb  —  .  —  —  r-  —  — —  -  ~  (3.25) 

•4m‘ 

For  two-dimensional  transverse  oscillations,  the  normal  mode  shape  ipm  is  a 
function  of  y  only.  Thus  Eq.  (3.2)  reduces  to 

(Plpri 


dy 2 


+  kmipm  =  0 


(3.26) 
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The  solution  for  the  fith  baffle  compartment  becomes 


0m  =  Acosk^y  +  £sink%ly 


(3.27) 


The  assumption  of  rigid  surfaces  gives  the  following  conditions  for  baffle  blades  and 
t’  chamber  wall  . 


d0m 

dy 


=  0 


y  n 
y  n 


H  =  1,2, 


(3.28) 


where  N  denotes  the  number  of  baffle  blades.  Thus,  the  normal  transverse  mode 
shape,  0m,  for  the  pth  baffle  compartment  is  determined  with  the  boundary 
conditions  set  by  Eq.  (3.28),  such  that 


0m  (y)  =  A  cos 


(3.29) 


If  the  baffle  blades  are  assumed  to  be  equally  spaced  in  the  chamber  and  have  rigid 
surfaces,  then  the  mode  shape  of  the  acoustic  pressure  in  the  /ith  baffle  compartment 
is  shown  to  be 


el?mMr  +Cmfce,9mt-2r 
e’9m»,iLk  -f 


(3.30) 


The  constant  denominator  within  the  square  bracket  is  added  in  this  expression 
to  facilitate  the  matching  at  the  interface.  This  constant  denominator  does  not 
affect  the  solution,  since  the  amplitude  of  linear  oscillations  is  arbitrary  anyway. 
An  extension  of  the  above  expression  to  cases  involving  a  non-uniform  distribution 
of  baffle  blades  is  straightforward,  only  at  the  expense  of  more  mathematical 
manipulations. 
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3.2.2  Acoustic  Field  in  Main  Combustion  Chamber 


The  solution  to  Eq.  (3.21)  for  the  axial  part  of  the  acoustic  field  in  the  main 
chamber  takes  the  form 


C„.(x)  = 


where 


MQ 


MQ 


MQ 


a 


t2 


+  6 2 

'  nc 


Qnc,1  -z 

a 


MQ 


a 


t2 


+  e 


(3.31) 


(3.32a) 


(3.326) 


Here,  the  subscript  c  denotes  quantities  associated  with  the  main  chamber.  The 
acoustic  boundary  condition  provided  by  the  exhaust  nozzle  may  be  represented  by 
a  nozzle  admittance  defined  below.  Although  this  is  a  simplistic  representation  of 
the  effect  of  the  nozzle,  the  ease  of  application  allows  for  a  focus  of  attention  on 
the  effect  of  baffles.  Furthermore,  for  short  nozzles,  this  representation  can  be  quite 
accurate.52 

t V  In  1 

(3.33) 


A  = 

N  p'/irp 


x=L 


Cnc  = 


(3.34) 


This  nozzle  admittance  function  leads  to  the  following  relationship  for  the  coeffi¬ 
cients  Anc  and  Bnc. 

_  B*C  _  9„.,1  +  ttpAj  _  ^  +  *».*  + 

’  «-«+n^  rn  _  +  J.c»  +  aA„ 

The  transverse  mode  function  xj)n  for  the  main  chamber,  based  on  the 
assumption  of  a  rigid  chamber  wall,  can  be  expressed  as 

^n(y)  =  Acos(^y) 


(3.35) 
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Substitution  of  Eqs.  (3.31)  and  (3.35)  into  Eq.  (3.1)  yields  the  expression  for  the 
mode  shape  of  the  acoustic  pressure  in  the  main  chamber. 

e‘9nc.,(*-i)  +  C„cei9"‘.’(x-L) 


=  cos(^)  4^ 

n=0  x  'I 


ALb-L) 


+  Cnce 


^^.ALb-L) 


(3.36) 


Again  the  denominator  in  the  big  bracket  is  added  for  the  convenience  of  matching 
at  the  interface  in  the  following  section. 

3.2.3  Matching  of  Acoustic  Fields  in  Baffle  Compartments 
and  Main  Chamber 


Together,  Eqs.  (3.30)  and  (3.36)  represent  the  acoustic  pressure  fields  in  the 
whole  chamber  in  terms  of  the  transverse  eigenfunctions.  Once  the  coefficients 
and  B„  and  frequency  fi  are  known,  the  pressure  distribution  in  the  entire  chamber 
can  be  easily  determined.  The  next  step,  therefore,  is  to  get  these  coefficients,  using 
the  matching  of  the  acoustic  fields  in  the  baffle  compartments  with  that  of  the  main 
chamber  at  the  interface.  This  matching  requires  continuity  of  acoustic  pressure 
and  axial  velocity. 


u 


x=L(,  P  1  x=Lb 

(3.37) 

,  =  ti'c|  , 

(3.38) 

From  Eqs.  (3.30)  and  (3.36),  the  pressure  distributions  of  the  baffle  compartments 

and  the  main  chamber  at  the  interface  (i.e.  x  =  Lb)  axe,  respectively 

-it  ^  (mzN  \ 

p  =  2^  Amcos[-^-y)^ 


\  h 

pc =s>sc°s(:fi') 

n=0  v  ' 


m= 0 
oo 


(3.39) 


(3.40) 


Note  that  indices  m  and  n  Eire  used  for  the  acoustic  fields  in  the  baffle  compartments 
and  the  main  chamber,  respectively,  to  avoid  confusion.  Substitution  of  the  above 
two  expressions  into  the  pressure  matching  condition  (3.37),  gives  an  expression  for 
the  baffle  compartment  coefficients,  A in  terms  of  the  main  chamber  coefficients, 


=o  r  cos(if y)d« 

<  =  £  - 

"=0L  ^(=Ty)dy 


(3.41) 


With  the  aid  of  the  acoustic  momentum  equation,  the  second  matching  condition 


(3.42) 


(3.38)  can  be  simplified  to  the  continuity  of  the  gradient  of  acoustic  pressures. 

dp ^  dpc 

it  =it  (3-42) 

OX  z=Lh  OX  z=Lh 

Now  taking  the  spatial  derivatives  of  (3.30)  and  (3.36)  with  respect  to  x,  and 
evaluating  ti  e  results  at  x  —  Lb  yields 

dp*  I  _ 


=  ^mCOsI  JT~y 

r=U  m^o  .  \  H  J 


(3.43) 


—  Lb  n=0  N  7 


eVdb-t)  +  Cnce,,"‘.’(L‘_L) 


(3.44) 


Substitution  of  Eqs.  (3.43)  and  (3.44)  into  Eq.  (3.42)  yields  an  explicit  expression 


for  the  coefficients  B„  in  terms  of  A independent  of  Eq.  (3.41). 

~  at  r  sT-uh  c™{!S7ry)cos{’fiy)dy 

Rc  _  V  V  _  x  -W1  “i 

n~hh[  m  LHcosH^y)dy  Hnm 


(3.45) 
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where 


Gm(fi)  = 


Hn(Q)  = 


l(ln 


+  iqm>.2Cm  beiq"".*Lb 


iie<c.l(^-L)  +  ^2Cnce.?n0(U- 


X) 


(3.46a) 


(3.466) 


eiqncA(U-L)  +  Cnce,q^iLi~L) 

Since  the  amplitude  is  arbitrary  for  linear  analyses,  a  normalization  of  the  acoustic 
pressure  to  a  particular  mode  within  the  main  chamber  helps  to  determine 
the  complex  frequency  of  acoustic  oscillations.  This  condition  can  be  formally 
expressed,  from  Eq.  (3.45),  as 


c»  N 

urn  =  £  £ 

m= 0  fi=l  L 


/<  1)H  cos(¥y)dv 

At  ^  - - — - x  Gm(Q) 


(3.47) 


f0  cos2(^y)dy 

where  h  is  the  specific  mode  of  concern  in  the  main  chamber.  For  example,  if  a  first 
transverse  mode  is  of  interest,  n  would  equal  one.  Note  that  this  normalization  is 
not  a  necessary  step  in  solving  the  problem,  but  it  is  only  a  convenience  used  to  aid 
in  the  numerical  computation  of  the  complex  frequency  and  mode  shapes.  If  this 
normalization  is  not  specified,  inversion  of  a  rather  large  matrix  would  be  required. 
This  matrix  contains  the  coefficients  of  Eqs.  (3.41)  and  (3.45),  with  eigenvalues 
corresponding  to  the  complex  frequency,  as  shown  in  the  Appendix.  Arbitrarily 
setting  the  dominant  mode  amplitude  equal  to  one  allows  a  simpler  iterative  solution 
to  the  system  of  equations,  and  will  not  change  the  value  of  frequency  or  mode 
shapes,  since  the  arbitrary  amplitude  cancels  out,  anyway. 

Because  of  the  transcendental  nature  of  Eq.  (3.47),  an  iteration  scheme  is 
used  to  determine  the  complex  frequency.  The  overall  numerical  procedure  is 
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conducted  as  follows.  First,  Eq.  (3.41)  is  substituted  into  Eqs.  (3.45)  and  (3.47)  for 
simplification, 
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The  second  step  involves  the  application  of  the  Newton- Raphson  technique  to 
determine  the  complex  frequency  from  Eq.  (3.49).  As  a  first  approximation, 
Bcm,  —  6m<n  and  Cl  =  u>„  are  used  to  initiate  the  iteration  procedure.  The 
calculated  Cl  is  used  to  update  B„  from  Eq.  (3.48),  which  is  then  substituted  back 
into  Eq.  (3.49)  to  calculate  the  complex  frequency.  This  procedure  repeats  itself 
until  the  converged  solutions  for  Cl  and  B„  are  obtained.  Finally,  the  coefficient 
is  determined  from  Eq.  (3.41).  The  iteration  scheme  functions  quite  effectively  and 
provides  reasonable  results  within  only  a  few  iterations. 
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3.3  Three-Dimensional  Baffled  Chamber 

In  this  section,  a  three-dimensioneil  cylindrical  baffled  chamber,  which  is  more 
representative  of  actual  rocket  systems,  is  investigated.  Usually  two  types  of  baffle 
configurations  are  used  to  suppress  the  transverse  modes  of  oscillation  in  such 
chambers:  radial  baffle  blades  and  hub-type  blades.  The  radial  baffle  blade  is 
utilized  to  control  tangential  modes  of  oscillation,  while  the  hub-type  affects  the 
radial  modes.  But  in  order  to  concentrate  attention  on  the  effect  of  baffles  on 
the  tangential  mode  of  oscillation,  only  problems  involving  simple  configurations  of 
radial  baffles  are  considered  here.  The  acoustic  field  can  be  treated  following  the 
general  approach  described  in  the  preceding  section. 


3.3.1  Acoustic  Field  in  Baffle  Compartments 

The  solution  of  Eq.  (3.21)  for  the  three-dimensional  chambers  takes  the  form 

Cm6(z)  =  A.mbetq”''> •» 1  +  Bmbetq”"‘-'x  (3.50) 

where  qmhA  and  qmi  3  are  defined  in  the  same  way  as  in  Eqs.  (3.23a)  and  (3.23 b). 
However,  the  transverse  mode  shape  ipn  for  a  cylindrical  chamber  with  radial  baffle 
blades  is  a  function  of  r  and  6  only,  taking  the  form 


ipn  =  (A  cos  m6  +  B  sin m0)Jm(Kmsr) 


(3.51) 


Baffle  blades  are  assumed  to  be  rigid  surfaces,  which  gives  the  following  conditions 
in  the  pth  baffle  compartment. 
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where  N  is  the  number  of  radial  baffle  blades.  These  boundary  conditions  are 
applied  to  Eq.  (3.51)  to  yield 


B  =  0 
nN 

m  =  —  (n  =  0,1,2,  •••) 


(3.53) 


In  most  actual  cases,  a  baffle  configuration  has  an  odd  number  of  radial  baffle  blades 
with  the  exception  of  a  one-blade  case,  which  is  practically  never  observed.  It  would 
offer  protection  from  modes  that  are  of  an  order  less  than  the  number  of  blades  and, 
to  some  degree,  protection  from  modes  of  an  order  higher  than  the  number  of  blades 
if  the  order  of  the  mode  divided  by  the  number  of  blades  is  not  equal  to  an  integer. 
With  the  odd  number  N,  then, 

Since  half-mode  waves  do  not  exist  physically,  the  only  values  of  m  considered  in 
this  analysis  are  those  for  which 


m  =  Nl  (l  =  0, 1,2,  •  •  •) 


(3.54) 


Thus,  the  transverse  normal  mode  shape  for  the  fit h  baffle  compartment  becomes 

<„=Acos(itfW)J„,(K„,.r)  (3.55) 

Here,  ipn  is  replaced  by  in  Eq.  (3.55)  to  distinguish  between  the  main  chamber 
and  the  baffle  compartment  indices. 

Substitution  of  Eqs.  (3.51)  and  (3.55)  into  Eq.  (3.1)  yields  the  mode  shape  of 
the  acoustic  pressure  for  the  /xth  baffle  compartment. 
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The  coefficient  Cmb  represents  the  boundary  condition  at  the  injector  face  and  is 
defined  in  the  same  way  as  in  Eq.  (3.25).  The  denominator  in  the  bracket  is  again 
added  for  convenience  of  matching  in  the  following  section. 


3.3.2  Acoustic  Field  in  Main  Combustion  Chamber 

The  solution  to  Eq.  (3.21)  for  the  acoustic  field  within  the  main  chamber  takes 
the  form 

Cnc(x)  =  4-  Bnct'q^*-L)  (3.57) 

where  qne  l  and  qnc  i  are  defined  in  the  same  way  as  in  Eqs.  (3.32a)  and  (3.326). 
The  normal  mode  shape,  t/>n,  for  standing  transverse  motion  in  the  main  chamber 
takes  the  form 

~  A  cos  ndJn(Kntr)  (3.58) 


To  avoid  confusion  of  the  main  chamber  solution  with  the  solution  in  the  baffle 
compartments,  the  index  nt  is  used  for  the  acoustic  field  in  the  main  chamber. 
Now,  combine  Eqs.  (3.57)  and  (3.58)  to  obtain  the  mode  shape  of  the  acoustic 
pressure  in  the  main  chamber. 
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'  cw«e,,(*-i)  +  Cnccif-c*(*"t) 
4-  Cnce’g"'-i(Lb~L) 


(3.59) 


The  coefficient  Cnc  incorporates  the  effect  of  the  nozzle  admittance  boundary 
condition  and  is  determined  in  the  same  way  as  in  Eq.  (3.34).  Again  the 
denominator  in  the  bracket  is  added  to  this  expression  to  facilitate  the  matching  in 
the  following  section. 
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3.3.3  Matching  of  Acoustic  Fields  in  Baffle  Compartments 
and  Main  Chamber 


Together,  Eqs.  (3.56)  and  (3.59)  express  the  acoustic  pressure  distributions 
in  the  whole  chamber  in  terms  of  eigenfunctions.  Once  the  coefficients  and 
B^t  the  frequency  are  found,  the  acoustic  fields  in  the  entire  chamber  can  be 
easily  computed.  The  matching  of  the  acoustic  fields  in  the  baffle  compartments 
and  the  main  chamber  at  the  interface  gives  sufficient  equations  for  determining 
the  coefficients  and,  thus,  the  wave  characteristics  of  the  baffled  chambers.  The 
matching  requires  continuity  of  acoustic  pressure  and  axial  velocity. 
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From  Eqs.  (3.56)  and  (3.59),  the  pressure  distributions  of  the  baffle  compartments 
and  the  main  chamber  at  the  interface  (i.e.  x  =  L&)  are,  respectively 
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Substitution  of  Eqs.  (3.62)  and  (3.63)  into  the  pressure  matching  condition  (3.60), 
yields  an  explicit  expression  for  the  baffle  compartment  coefficients,  A^u ,  in  terms 
of  the  main  chamber  coefficients,  B„t. 
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With  the  aid  of  the  acoustic  momentum  equation,  the  second  matching  condition, 
Eq.  (3.61),  can  be  interpreted  as  the  continuity  of  the  axial  gradient  of  acoustic 


pressure. 
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(3.65) 


Differentiation  of  Eqs.  (3.56)  and  (3.59)  with  respect  to  x  and  subsequent  evaluation 
at  x  =  Lt,  give 
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Substitution  of  Eqs.  (3.66)  and  (3.67)  into  Eq.  (3.65)  gives  an  expression  for  the 
Fourier  coefficients  of  the  main  chamber,  BcnV  in  terms  of  those  coefficients  for  the 
baffle  compartments,  . 


oo  oo  N 
1=0  3=1  *1=1  L 


2*e.  R 

f  I  V’v,  Pnti'drdd 

J  **■(*-!)  Jt>  j-, 

ah  n _ ~  ^jyi 

NU  r3w  rR  t  ,  .  \2 


Xty 


//  J0  {PntYrdrdO  Entity  J 


(3.68) 


where 


G„,.(ty  = 


Entity  = 


*9n 


iq 


nc,  1 


+  Cmhtxq”'i’3Lb 

+  iqne,Cnctiq^Lb-L) 
+  Cnce*,nc-j(L*~L) 


(3.69a) 

(3.695) 


As  was  done  for  the  two-dimensional  case,  a  normalization  of  the  acoustic  pressure 
to  a  particular  mode  within  the  main  chamber  is  made  in  determining  the  complex 
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frequency  of  the  acoustic  oscillations.  This  normalization  can  be  expressed,  from 
Eq.  (3.68),  as 
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where  n  and  t  are  specific  tangential  and  radial  modes  of  concern,  respectively,  in 
the  main  chamber.  For  example,  if  a  first  tangential  mode  is  under  consideration, 
both  n  and  t  would  equal  one. 

An  iteration  scheme  is  again  used  to  determine  the  complex  frequency,  and 
the  numerical  iteration  proceeds  as  follows.  First,  substitution  of  Eq.  (3.64)  into 
Eqs.  (3.68)  and  (3.70)  gives 
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Second,  a  Newton-Raphson  iteration  scheme  is  applied  to  determine  the  complex 
frequency  ft  from  Eq.  (3.72).  As  a  first  approximation,  and  ft  =  u>„ 
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axe  used  to  initiate  the  iteration  procedure.  The  calculated  Q  is  used  to  update 
B„t  from  Eq.  (3.71),  which  is  then  substituted  back  into  Eq.  (3.72)  to  compute  the 
complex  frequency.  This  procedure  is  repeated  until  converged  solutions  for  Q,  and 
B„t  are  obtained.  Finally,  the  coefficients,  A*  are  determined  from  Eq.  (3.64). 
Like  the  two-dimensional  solution,  this  method  gives  reasonable  values  for  complex 
frequency  Oina  few  iterations. 

3.4  Results  and  Discussion  of  Linear  Acoustic  Analysis 

The  analysis  described  in  the  previous  sections  has  been  applied  to  the  study 
of  acoustic  wave  characteristics  in  baffled  combustion  chambers,  and  the  results 
have  provided  useful  information  about  the  effect  of  baffles  on  the  detailed  flow 
structure.  The  consequences  for  the  two  cases  are  presented  separately:  two-  and 
three-dimensional  baffled  chambers. 

3.4.1  Acoustic  Field  in  Two-Dimensional  Baffled  Chamber 

To  gain  confidence  in  the  accuracy  of  this  analysis,  calculated  acoustic  fields 
of  normal  modes  axe  first  validated  with  numerical  solutions  obtained  from  a 
well-established  finite  element  code,  ANSYS.53  The  ANSYS  code  uses  a  3-D  iso¬ 
parametric  acoustic  fluid  element  and  has  been  successfully  used  to  study  sound 
propagation  and  submerged  structure  dynamics. 

Figure  3.3  shows  the  acoustic  pressure  fields  of  the  first  transverse  mode  in 
chambers  containing  one  and  two  baffle  blades.  The  aspect  ratio  (L/H)  of  the 
chamber  is  1.25,  similar  to  that  of  the  two-dimensional  ‘equivalent’  combustion 
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Acoustic  Pressure  Contours  in  Baffled  Chambers: 
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chamber  of  the  F-l  engine.  Zero  acoustic  admittance  is  specified  at  the  injector 
and  nozzle  entrance  to  ensure  normal  mode  solutions.  Because  of  the  limitation 
of  the  wave-front  number  in  the  ANSYS  code,  small  kinks  can  be  observed  in  the 
numerically  calculated  pressure  contours.  The  analytical  solutions  are  in  excellent 
agreement  with  these  finite  element  calculations,  so  this  analytical  approach  can  be 
applied  to  acoustic  problems  in  the  baffled  combustion  chamber  with  confidence. 

Figure  3.4  presents  the  effect  of  baffle  length  on  the  oscillatory  pressure  fields  of 
the  first  transverse  mode  in  a  one-baffle  chamber,  and  shows  that  the  baffle  exerts 
significant  influence  on  the  acoustic  field.  For  long  baffles  (e.g.,  Lb  =  0.4/f ),  the 
acoustic  motion  inside  the  baffle  compartments  very  much  resembles  a  longitudinal 
wave  and  transits  to  a  transverse  oscillation  downstream  of  the  baffle.  This  is  seen 
clearly  by  noticing  that  the  pressure  nodal  line  for  the  first  transverse  mode  in  the 
unbaffled  chamber  lies  at  the  center  of  the  chamber  parallel  to  the  chamber  wall.  But 
once  the  baffle  is  added  to  the  chamber,  it  must  move  to  satisfy  the  acoustic  velocity 
boundary  conditions  set  by  the  baffle  blades,  such  that  the  pressure  contour  lines 
are  normal  to  the  baffle  surface.  By  shifting  the  mode  of  oscillation  from  transverse 
to  longitudinal  inside  the  baffle  compartments,  baffles  can  decrease  the  prevalence 
of  the  transverse  modes,  which  have  historically  been  most  troublesome  in  liquid 
rocket  combustion  instability  problems,  whereas  longitudinal  modes  have  not.  As 
a  result,  however,  regions  of  high  amplitude  pressure  fluctuation  are  concentrated 
near  the  injector  face.  If  pressure  sensitive  combustion  processes  are  responsible  for 
the  instability,  the  addition  of  a  baffle  might,  therefore,  worsen  the  problem,  since 
the  region  sensitive  to  the  variation  of  acoustic  motion  is  usually  close  to  the  injector. 
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The  concentration  of  large  amplitude  pressure  fluctuations  in  this  region  thus  has 
the  potential  to  release  large  amounts  of  energy.  However,  in  an  unbaffled  chamber, 
the  pressure  variation  associated  with  a  transverse  mode  can  change  the  orientation 
of  the  propellant  jets  issuing  from  the  injector.  This  can  have  a  significant  influence 
on  the  atomization  and  spray  characteristics  near  the  injector  face  that  may  be 
coupled  with  the  chamber  acoustics  to  generate  instability.  With  the  addition  of 
a  baffle,  however,  the  longitudinal  nature  of  pressure  oscillations  near  the  injector 
face  might  have  no  effect,  or  at  least  a  substantially  different  effect  on  this  type  of 
combustion  process  near  the  injector  face. 

Figure  3.5  shows  the  acoustic  pressure  contours  of  the  first  transverse  mode 
in  a  two-baffled  chamber.  For  a  short  baffle  length  (e.g.,  L\,  =  0.1  H),  transverse 
oscillations  exist  in  almost  the  entire  chamber  except  for  regions  near  the  baffle  tips, 
but  the  degree  of  wave  distortion  increases  with  increased  baffle  length.  As  with  the 
one-baffle  case,  longitudinal  oscillations  become  conspicuous  in  the  upper  and  lower 
baffle  compartments  for  X*  >  0.2 H,  which  might  mean  enhanced  stability,  as  argued 
previously.  Note  that  the  acoustic  pressure  is  almost  uniformly  distributed  in  the 
center  baffle  compartment,  with  its  magnitude  at  a  minimal  level.  This  suggests 
that  the  center  compartment  may  be  the  most  effective  in  separating  the  acoustic 
waves  from  sensitive  combustion  processes,  and  thus  might  be  able  to  overcome  any 
possible  destabilizing  influences  from  the  adjacent  compartments. 

However,  the  situation  for  the  second  transverse  mode  is  quite  different.  Figure 
3.6  shows  the  oscillatory  pressure  fields,  which  indicates  that  the  wave  motion 
exhibits  longitudinal  motion  in  the  top  and  bottom  baffle  compartments,  similar 
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to  the  first  mode  case.  However,  large-amplitude  bulk  oscillations  now  exist  in  the 
center  region.  If  pressure  sensitive  processes  are  responsible  for  the  instability, 
this  concentrated  pressure  fluctuation  could  promote  instability.  The  effect  of 
acoustic  velocity  on  unsteady  combustion,  however,  is  quite  limited  in  the  baffle 
compartments,  especially  in  the  center  one.  According  to  the  experiments  of  Levine 
and  Bambanek38  and  Reardon,35  transverse  components  of  the  acoustic  velocity 
have  a  strong  effect  on  the  combustion  process.  For  instance,  atomization,  local 
mixture  ratio,  and  adjacent  injector  element  interactions  can  all  be  affected  by 
acoustic  velocity  fluctuations  in  the  transverse  direction.  Near  the  injector  face, 
where  large  temperature  and  concentration  gradients  are  found,  this  effect  cam 
be  especially  pronounced.  From  Fig.  3.7,  however,  it  is  clear  that  the  transverse 
components  are  very  small  near  the  injector  face  inside  the  baffle  compartments. 
Figure  3.7  shows  the  acoustic  pressure  and  velocity  vector  together  for  the  first 
transverse  mode  in  a  two-baffle  chamber. 

So  for  an  unbaffled  chamber  in  which  a  transverse  mode  of  instability  arises 
from  velocity  sensitive  processes,  addition  of  baffles  will  likely  have  a  stabilizing 
effect.  Thus,  the  longitudinalization  of  transverse  oscillations  brought  about  by 
the  baffle  can  be  supposed  to  have  stabilizing  or  destabilizing  effects,  depending 
primarily  on  the  relative  importance  of  pressure  sensitive  and  velocity  sensitive 
combustion  processes. 

The  second  transverse  mode  for  the  one-baffle  chamber  and  the  third  transverse 
mode  for  the  two-baffle  chamber  are  also  examined,  as  shown  in  Fig.  3.8.  The 
acoustic  fields  are  almost  identical  to  those  in  an  unbaffled  chamber,  showing 


Figure  3.7  Acoustic  Pressure  and  Velocity  for  the  First 
Transverse  Mode  in  a  Two-Baffle  Chamber 
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that  the  baffles  have  little  effect  on  these  transverse  modes  because  the  boundary 
conditions  imposed  by  the  baffle  blades  are  automatically  satisfied  if  they  are  located 
at  the  acoustic  pressure  anti-nodal  (or  the  velocity  nodal)  lines. 

The  effect  of  baffle  length  on  the  oscillation  frequency  is  also  studied. 
Comparison  with  the  numerical  ( ANSYS)  results  is  made  again,  since  no  appropriate 
experimental  data  was  found.  Figure  3.9  presents  the  results  for  the  first  transverse 
mode  in  a  one- baffle  chamber,  where  the  frequency  is  normalized  with  respect  to 
the  natural  frequency  of  an  unbaffled  chamber.  Excellent  agreement  is  observed 
between  the  analytical  predictions  and  the  finite-element  calculations,  with  relative 
errors  of  less  than  one  percent.  The  frequency  decreases  persistently  with  increased 
baffle  length,  perhaps  due  to  an  increase  in  the  effective  path  length  of  the  gas 
particles,  as  caused  by  the  baffle  blade.  This  same  trend  can  also  be  observed  in 
actual  three-dimensional  engine  tests.  For  a  baffle  length  of  0.4# ,  the  reduction 
of  frequency  from  that  of  the  unbaffled  chamber  is  greater  than  25%.  The  ensuing 
influence  on  the  liquid-propellant  combustion  response  can  be  significant  and  must 
be  taken  into  account  in  a  realistic  treatment  of  combustion  instabilities  in  baffled 
combustors.  Since  combustion  instability  depends  on  a  specific  relationship  between 
the  acoustic  frequency  and  the  characteristic  response  time  of  combustion  processes, 
sufficient  modification  of  the  oscillation  frequency  will  probably  lessen  the  amount 
of  acoustic/combustion  coupling,  and  might  terminate  the  instability  altogether. 
Thus  decreasing  oscillation  frequency  is  another  potential  stabilizing  effect  of  the 
baffles  on  the  chamber. 

It  should  be  noted  that  numerical  calculations  using  ANSYS  takes  several  hours 
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Figure  3.9  Effect  of  Baffle  Length  on  the  Frequency  of  the 
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on  the  VAX-8700  system,  but  only  a  few  minutes  is  required  to  get  analytical  results 
with  40  terms  in  the  eigenfunction  expansion.  This  suggests  that  this  analysis  has 
many  strong  features. 

As  mentioned  before,  a  critical  issue  in  the  study  of  liquid  rocket  combustion 
instability  is  the  effect  of  acoustic  velocity  on  propellant  combustion.  In  considera¬ 
tion  of  this  point,  the  acoustic  velocity  field  is  calculated  in  detail,  giving  the  result 
shown  in  Fig.  3.10.  The  acoustic  velocity  reaches  its  maximum  amplitude  due  to 
the  large  pressure  gradient  near  the  baffle  tip,  then  decreases  in  the  far  field.  The 
transverse  components  of  the  acoustic  velocity  are  generally  much  larger  than  the 
axial  components  just  downstream  of  the  baffle  tip,  since  the  acoustic  wave  front 
must  turn  180  degrees  around  the  baffle.  Figure  3.10(b)  shows  the  distributions  of 
transverse  acoustic  velocity  at  various  axial  locations,  where  (v')*  stands  for  the 
velocity  normalized  with  respect  to  its  maximum  value  in  the  transverse  plane  at 
the  given  axial  location.  In  the  baffle  compartment  (z  <  0.4/T),  the  requirement 
for  transverse  acoustic  velocity  to  vanish  at  the  baffle  blade  causes  two  peaks  in 
the  velocity  profile.  The  profile  then  becomes  a  simpler  distribution  downstream 
of  the  baffle,  with  a  single  peak  at  the  centerline.  Note  that  the  magnitude  of  the 
transverse  velocity  components  near  the  injector  face  is  negligibly  small  due  to  the 
longitudinalization  of  the  wave  there. 

Near  the  baffle  tip,  however,  large  magnitudes  and  a  rapid  variation  in  the 
velocity  profile  are  observed.  Furthermore,  the  large  gradients  of  acoustic  pressure 
observed  near  the  baffle  tip  provide  sizable  centripetal  forces  as  the  wave  turns 
180  degrees  around  the  baffle  tip  in  satisfaction  of  the  geometric  condition  set  by 
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the  baffle.  Baer  and  Mitchell31  claimed  that  viscous  dissipation  of  acoustic  energy 
due  to  these  rapid  variations  of  acoustic  velocity  is  crucial  in  stabilizing  the  baffled 
rocket  engine.  It  is  worth  noting  here  that  the  drastic  change  of  acoustic  velocity 
near  the  baffle  tip  should  lead  to  the  formation  of  vortices,  which  are  convected 
downstream  in  the  presence  of  mean  flow.  In  the  case  of  oscillatory  flow  past 
a  baffle  compartment,  it  might  be  expected  that  vortices  are  also  created  in  the 
compartments,  with  the  direction  of  rotation  alternating  in  response  to  the  flow 
past  the  baffle.25  Such  creation  and  destruction  of  vortices  may  provide  for  energy 
dissipation  that  is  significant  in  certain  cases.  However,  the  longitudinalization  of 
flow  in  the  compartments  suggests  that  downstream  shedding  of  vortices  is  more 
important,  especially  in  the  case  of  the  longer  baffles.  The  subsequent  influence  of 
these  vortices  on  the  mean  flowfield  may  also  significantly  alter  the  combustion 
distribution  and  instability  characteristics,  and  is  addressed  by  Wicker32  in  a 
systematic  fashion. 

3.4.2  Acoustic  Field  in  Three-Dimensional  Baffled  Chamber 

Calculations  are  also  carried  out  for  acoustic  fields  in  three-dimensional 
cylindrical  chambers  with  radial  baffles,  thereby  simulating  more  realistic  rocket 
engine  configurations. 

Figure  3.11  presents  the  normalized  frequency  of  the  first  tangential  mode 
as  a  function  of  the  baffle  length  and  the  number  of  blades.  Good  agreement 
with  experimental  data20'22,54,55  was  obtained.  As  with  the  two-dimensional  case, 
the  oscillation  frequency  decreases  with  baffle  length  due  to  the  increased  acoustic 
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Figure  3.11  Effect  of  Baffle  Length  on  the  Frequency  of  the 
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path  length  around  the  baffle  blades.  This  depression  of  the  oscillation  frequencies 
can  eliminate  acoustic/combustion  coupling  in  the  chamber,  or  perhaps  alter  the 
chamber  mode  so  that  the  instability  may  be  avoided.  This  phenomenon  might 
serve  as  one  of  the  stabilizing  effects  of  the  baffles.  It  is  also  interesting  to  mention 
that  the  number  of  baffle  blades  plays  a  minor  role  in  determining  the  frequency. 
Increasing  the  number  of  blades  from  three  to  seven  only  slightly  decreases  the 
frequency,  a  phenomenon  also  supported  by  experimental  data.  This  aspect  is 
understandable  considering  that  when  more  than  three  baffle  blades  are  used,  the 
first  tangential  mode  is  always  affected,  regardless  of  the  number  of  blades. 

Figure  3.12  shows  the  acoustic  pressure  contours  at  various  cross  sections  in 
a  chamber  with  three  radial  baffle  blades.  This  baffle  configuration  is  the  simplest 
one  that  can  get  rid  of  first  tangential  mode  instability,  which  is  the  most  often 
seen  and  the  most  troublesome  in  the  actual  rocket  systems.  The  wave  resembles  a 
classical  tangential  mode  in  the  downstream  end  of  the  chamber  (i.e.  x/L  =  1.0), 
but  the  pressure  contours  become  distorted  near  the  baffle  tips,  and  finally  reach 
a  symmetric  pattern  dictated  by  the  geometric  constraints  of  the  baffle  blades.  As 
seen  for  the  two-dimensional  case,  unsteady  motions  in  the  injector  end  of  the  baffle 
compartments  have  a  more  longitudinal  or  even  bulk  type  oscillation,  with  a  quite 
uniform  transverse  spatial  distribution. 

Also  important  is  that  the  pressure  amplitude  is  higher  than  that  in  the  main 
chamber.  This  peculiar  phenomenon  suggests  that  the  baffle  tends  to  redistribute 
the  acoustic  energy  in  the  chamber  from  the  downstream  to  the  upstream  region. 
Usually  most  of  the  energy  release  from  liquid  propellant  occurs  in  the  upstream 
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region.  The  high-amplitude  acoustic  pressure  in  that  region  has  the  potential  to 
release  larger  amounts  of  energy,  and  could  have  a  net  destabilizing  effect  on  the 
engine.  This  mechanism  is  formalized  in  terms  of  Rayleigh’s  criterion,29  which 
states  that  if  a  pulse  of  heat  is  added  at  a  time  and  location  where  the  acoustic 
pressure  is  increasing,  the  pressure  rise  is  encouraged  and  hence  the  heat  addition 
has  a  destabilizing  effect.  However,  as  was  already  noted  in  two-dimensional 
results,  the  characteristics  of  the  wave  motions  in  the  baffle  compartments  are 
significantly  changed.  Under  further  considerations,  this  pressure  concentration 
may  not  destabilize  the  engines. 

Figure  3.13  presents  the  transverse  acoustic  velocity  vectors  at  various  cross 
sections.  The  acoustic  velocities  in  the  baffle  compartments  have  very  small  ampli¬ 
tudes,  especially  near  the  injector  face,  where  large  temperature  and  concentration 
gradients  exist.  Since  the  transverse  components  of  acoustic  velocity  would  be 
expected  to  have  little  effect  on  the  combustion  processes  in  that  region,  addition 
of  baffle  blades  in  the  chamber  could  yield  a  stabilizing  effect,  especially  if  velocity 
sensitive  combustion  processes  are  of  primary  concern.  Large-amplitude  velocity 
oscillations  exist  immediately  downstream  of  the  baffles,  which  is  similar  to  the 
two-dimensional  case.  The  large  acoustic  velocity  gradients  in  this  region  lead 
to  the  formation  of  vortices,  which  are  convected  downstream  by  the  mean  flow. 
Viscous  dissipation  within  each  vortex  and  convection  of  energy  associated  with 
the  vortical  structure  out  of  the  chamber  may  be  a  significant  mechanism  of  energy 
dissipation.32  Also,  the  influences  of  these  vortices  on  the  mean  flowfield  might 
change  the  combustion  distribution  and  instability  characteristics. 
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CHAPTER  4 

NONLINEAR  SPINNING  TRANSVERSE  OSCILLATIONS 


Before  dealing  with  the  nonlinear  acoustic  analysis  in  the  baffled  combustion 
chamber,  this  chapter  addresses  the  characteristics  of  spinning  wave  motions  in 
a  cylindrical  unbaffled  chamber.  Spinning  transverse  oscillations  are  quite  often 
observed  in  many  types  of  combustion  chambers  and  are  known  to  be  detrimental 
due  to  the  increasing  combustion  efficiency  through  the  acceleration  of  mixing 
processes.  The  purpose  of  this  chapter  is  to  examine  how  the  spinning  wave  motions 
are  related  to  the  velocity  coupling,  and  the  conditions  for  which  stable  limit  cycles 
exist  for  the  special  cases  of  combinations  of  spinning  modes. 


71 


A  set  of  second-order  ordinary  differential  equations  for  the  time-varying 
amplitude  of  each  mode  can  be  obtained  by  averaging  the  conservation  equations 
over  the  volume  of  the  chamber  using  Galerkin’s  method.  Now  multiply  Eq.  (2.29) 
by  V’n  and  Eq.  (4.2a)  by  p' ,  subtract  the  results,  integrate  over  the  volume,  and 
apply  Green’s  theorem  to  find 

+  kl^nP')dV  =  -JJJ % PnhdV  -  <^Bf  •  ndS  (4.3) 

Substitution  of  Eq.  (4.1)  into  Eq.  (4.3)  and  some  rearrangements  of  the  results  lead 
to  a  system  of  second-order  ordinary  differential  equations  for  the  time-dependent 
amplitudes. 

+  u2nT}n  =  Fn  (4.4) 


where  the  forcing  function  Fn  can  be  expressed 

F"=w?\SlhMV+fJndS. 

and 


n2dV 


(4.5a) 


(4.56) 


The  next  step  is  to  represent  the  forcing  function  Fn  in  terms  of  time-varying 
amplitudes  In  this  chapter,  the  main  focus  is  on  the  influence  of  the  linear 

interactions  between  modes  on  the  spinning  wave  motions  within  the  chamber,  so 
that  only  the  linear  part  of  the  forcing  function  is  taken  into  consideration.  The 
linear  approximation  of  /13  in  Eq.  (2.30c)  is 


h 


3  tin 


R  A  Hv  d  . , 
d2  Cv  dtu‘ 
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where  AHV  is  the  heat  of  combustion  of  the  liquid  propellant  per  unit  mass.  This 
expression  is  not  complete  until  the  burning  rate  is  represented  in  terms  of  flow 
variables.  Such  closure  of  the  solution  is  conveniently  accomplished  through  the 
use  of  the  combustion  response  function  concept.  The  original  definition  of  the 
combustion  response  function  uses  only  pressure  for  representing  the  burning  rate, 
but  according  to  Reardon,35  the  transverse  velocity  components  play  an  important 
role  on  the  burning  rate  fluctuation  for  the  spinning  mode  of  oscillations.  In  this 
analysis,  the  velocity  coupling  effect,  as  well  as  the  pressure  coupling,  is  incorporated 
with  the  combustion  response  function  to  express  the  burning  rate.  The  unsteady 
combustion  is  modeled  as 


A,  _  p  _  v  _  w' 
—  —  77.p  —  +  7lv  —  7ZW  — 


wi 


(4.7) 


where  v'  and  w'  are  acoustic  velocity  perturbations  in  the  radial  and  tangential 
directions,  respectively;  the  pressure- coupled,  radial  and  tangential  velocity-coupled 
response  functions  are  symbolized  by  1ZP,  1lv  and  1ZW,  respectively.  For  the  linear 
acoustic  problems,  all  variables  are  assumed  to  vary  in  a  time-harmonic  fashion, 

eg-, 

p'(r,()  =  p(r)e'n‘ 

(4.8) 

Since  the  combustion  response  function  should  account  for  both  the  amplitude  and 
the  phase  relationship  between  flow  variables  and  burning  rate  perturbation,  the 
complex  expressions  for  7 Zp,  7ZV  and  7ZW  are  required,  thus,  Eq.  (4.7)  is  rewritten 
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as  the  sum  of  its  real  and  imaginary  parts. 


oj 

u>i 


+  Tle(Tlv)~  +  ^-lm(Tlv)%- 
a  ail  at 

_  .it/  1  _  .  dw1 

+  7Ze(7Zw)~  +  —-Im(7lw)  — 
a  ail  at 


(4.9) 


Then,  the  linear  approximation  of  the  combustion  response  related  source  term 
(h3lin)CR  takes  the  form, 

1 
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„  d2v' 
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d2wr 

1  dt 

+  hn 

(4.10a) 


where, 
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(4.106) 


hn  = 


R  u?i2m(7lw)AHv 


Cv  Ua 

Because  the  combustion  response  function  is  a  function  of  frequency,  all  coefficients 
in  Eq.  (4.106)  are  also  functions  of  the  frequency  and  have  different  values  for  each 
mode  of  oscillations. 
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If  the  series  expansions  for  u'g  and  p'  in  Eq.  (4.1)  and  the  linear  approximation 
of  combustion  response  related  source  term  (h3lin)CR  in  Eq.  (4.10a)  are  substituted 
into  Eq.  (4.5),  the  linear  part  of  the  forcing  function  then  becomes 

OO 

{F n)ltnear  —  ^  ^{DniPi  +  -Ent^i)  (4.11) 

i=l 

where  the  linear  coefficients  _Dni  and  Ent  axe  given  by, 


-  —  x  (V  x  uf )]  •  V</>„  +  yrpitpnV  •  us 


dxfri  ,  ,  I2n  1  dxfri  , 

Gin^i^n  +  — +  — - SZ’V’n 

p  or  p  r  Off 


dV  (4.12a) 


r  _  1  m\„  1  2  ;  /  -2  ,  50,  ,  ,  ^ln  1  50,  , 

En‘  ~  £n2  JjJ  [G2nk^iVn<1  +  ^  dr  +  ^  r  de 


dF  (4.126) 


Here,  the  linear  interaction  effects  between  modes  represented  by  the  coefficients  in 
Eq.  (i  106)  are  incorporated  with  the  linear  parameters  in  Eq.  (4.12). 

For  the  standing  mode  of  oscillations,  each  tangential  motion  has  only  one  mode 
shape  in  the  azimuthal  direction.  Therefore,  the  linear  coupling  between  modes 
cannot  exist,  even  though  the  linear  contribution  of  the  combustion  response  is 
considered.  For  the  spinning  mode  of  oscillations,  however,  two  periodic  motions  in 
a  tangential  direction  are  possible.  That  is,  the  combination  of  0n=cos  m9Jm(nmtr) 
and  0,=  sin  m9Jm(Kmar)  is  possible  for  the  spinning  transverse  oscillations.  For 
the  first  three  modes,  the  wave  numbers  and  mode  shapes  are  given  by, 


75 


First  Tangential  Mode  (IT) 


k.\\R  =  1.8412;  -0i  =  cos0Ji(Knr),  fa  =  sin0Ji(/cur)  (4.13a) 


First  Radial  Mode  (1R) 

«oi  R  =  3.8317;  fa  =  Jo(Koir)  (4.136) 

Second  Tangential  Mode  (2T) 


k21R  =  3.0542;  fa  =  cos  2072(*2ir)>  =  sin20J2(/C2i^)  (4.13c) 


The  linear  coefficients  and  E’ni  can  be  computed  from  Eq.  (4.12)  for  the  first 
three  modes.  Unlike  Yang  and  Culick’s  standing  mode  analysis,47  the  linear  coupling 
coefficients,  Du,  Eu,  D41,  and  E41  appear  for  the  first  tangential  mode  and  D35, 
E35,  H53  and  £53  for  the  second  tangential  mode,  with  the  following  values. 


Du  =  D44  =  -Gn  +  ?2-Ax 

P 


D22  —  — G12  H — r-  A2 


D33  =  D55  =  —G 13  + 


H ■ 


23 


D\4  =  —D41  =  ~rA4 

P 

D35  =  —D53  —  ~~-A5 

P 


(4.14a) 


And 


E\\  =  E44  =  G2\w\  +  — r-Ai 

P 

E22  =  G22<jJ2  d - - ~A2 
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Ezi  —  Ess  =  G2sU}\  4 - ~j43 

P 

Ei  4  =  —Eu  =  -~A4 
p 

Ess  —  —Es3  =  -~As 
P 


(4.146) 


where 


Ai  — 


(i  -  i£v)[MKuR)i 


J  rJiiitiar)  jj0(Ki,r)  -  ^(ki***)] dr  (4.15a) 


A2  =  — 


As  = 


A,  = 


- —  —  - - ~r  f  rJ0(Koar)Ji(Kosr)dr  (4.156) 

R2[Jo(kosR)]  Jo 

- 7 — y1 - rj  /  rJ2(K2ar)  Ji(K2ar)  ~ ’M'^*')]  <M4-!5c) 

i1  - -^wnM^sR)}  Jo  J 

jf  l- '*M  *  (415e) 


(4.15rf) 


(4.15e) 


An  approximate  solution  technique  using  the  method  of  time-averaging  gives 
the  temporal  behavior  of  the  amplitude.  First  of  all,  time- varying  amplitudes  rjn{t ) 
are  assumed  to  have  the  form, 


rjn(t)  =  rn(f)sin[u>nf  +  <f>n(t)] 


(4.16) 


where  rn(<)  and  <f>n(t )  are  the  amplitude  and  phase  of  the  limit  cycle,  respectively. 
Using  the  linear  coupling  coefficients  in  Eq.  (4.14)  and  applying  the  time-averaging 
method,  the  equations  for  time- varying  amplitudes  due  to  the  combustion  response 
are  given  by, 
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First  Tangential  Mode  (IT) 


dr\ 

(-jj-)cR  =<*i4r4  cos(<?!>i  -<f>A)~  01ArA  sin(«b  -  <f>A) 
dr  4 

(-jj-)cR  =«4in  cos(^!  -  <j)A)  -t-  eA 4r,  sin(^i  -  <t>4 ) 

(ri  ~^^CR  =  -  ai4ri  sin^1  “  ^  ~  °ur4  cos(^i  -  <M 

(r4^-)cR  =a4in  sin(^!  -  <^4)  -  cos(<£j  -  <j>A ) 


Second  Tangential  Mode  (2T) 


dr$ 

(~jf)cR  =«35r5  cos(^>3  -  cf> 5)  -  ^35^5  sin(<63  -  <£5) 

(— ij~)cR  ~a53r3  COs(<^»3  —  <£5)  +  $537*3  sin(03  —  <£5) 
(r3^-)CR  =  ~  «35r5  Sin(^3  “  &>)  “  ^35^5  COs(<£3  -  <£5) 
(r5^-)CR  =a53r3  sin(</>3  -  <j>5)  -  $53r3  cos(<£3  -  <f>5) 


where, 


<*ni  =  -Dni/2  9ni  =  -Eni/2u>n 


(4.17 a) 
(4.176) 
(4.17c) 
(4.17c?) 

(4.18a) 

(4.186) 

(4.18c) 

(4.18d) 

(4.19) 


4.2  Nonlinear  Spinning  Oscillations 

Equations  for  time-varying  amplitudes  due  to  the  combustion  response  in 
Section  4.1  are  combined  with  those  in  Ref.  47  to  investigate  the  nonlinear 
spinning  transverse  oscillations.  For  simplicity,  only  the  first  three  modes  (IT, 
1R  and  2T)-the  most  commonly  observed  in  practical  systems-are  considered  here. 
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Four  simple  combinations  of  modes  are  examined  for  the  spinning  wave  motions: 
the  first  tangential  mode,  the  first  tangential/first  radial  modes,  and  the  first 
tangential/second  tangential  modes.  After  this,  the  combination  of  three  modes 
is  discussed.  The  conditions  for  which  the  limit  cycles  exist  are  studied  for  each 
case,  and  the  influence  of  linear  coupling  parameters  on  the  spinning  wave  motions 
is  also  investigated. 


4.2.1  First  Tangential  Mode 


The  first  tangential  mode  alone  can  exhibit  the  spinning  wave  motion  if  the 
linear  coupling  terms  are  taken  into  account.  Equations  of  time- varying  amplitudes 
and  phases  due  to  the  combustion  response  (Section  4.1)  are  added  to  those  of 
transverse  oscillations  in  Ref.  47. 
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where, 


X  =  <f>i  —  <f>4 


The  relations  between  linear  parameters  are  easily  found  as  follows, 

ot i  =  o4 ,  6}  =  84 

a  14  —  — «41»  014  =  —041 


(4.20) 

(4.21) 


(4.22) 


(4.23) 


(4.24) 
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Similar  to  the  standing  mode  analysis,47  the  limiting  values  of  amplitude  rn  and 
phase  <j>n  are  given  in  the  following  form, 

rno  =  constant  (4.25a) 

<j>nQ  =  vnt  +  in  (4.256) 

The  subscript  0  denotes  the  value  in  the  limit  cycle.  Comb  ining  Eqs.  (4.23)  and 
(4.25)  gives 


v\  =  1/4  (4.26a) 

X0  =  ix  -  U  (4-266) 

At  the  limit  cycle,  then,  Eqs.  (4.20)-(4.22)  become 

0  =  aino  +  a^no  cos X0  -  0i4r4o  sinX0  (4-27) 

0  =  a^no  +  O4irio  cos  Xo  +  @4irjo  sinXo  (4.28) 

0  =  —{Q\  —  ^4) 


^14r40  #41^10  v  «14^40  .  «41»'10  .  v 

— - cosAo  — - 1 - smXo  (4.29) 

no  no  r10  r40 

The  practical  assumption  in  a  reasonable  sense  is  taken  for  the  limiting  amplitudes. 

no  =  no  (4.30) 

After  some  manipulations  of  Eqs.  (4.27)-(4.29),  the  relation  between  frequency 
modulations  is  obtained  as, 

i\=it  +  -ny--  (n  =  1,2,3  ■•  •)  (4.31) 
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This  means  that  mode  1  and  mode  4  have  a  phase  difference  of  either  j  or  at 
the  limit  cycle.  On  the  other  hand,  acoustic  pressure  for  the  first  tangential  mode 
is 


p'(r,t)/p  =  r10  sin(u>i<  +  <f>  10)  cos 0J\(^\r) 

+  r40  sin(u4t  +  (f>40)  sin  6 J1(n1r) 

Equation  (4.32)  can  be  further  simplified  using  Eqs.  (4.30)  and  (4.31), 


(4.32) 


p'(r,t)/p  =  r10sin[(ui  +  iq )t  +  6  ±8]MK\r)  (4.33) 


Equation  (4.33)  shows  that  the  spatial-  and  time- dependent  parts  of  the  first 
tangential  mode  are  combined  together,  representing  the  spinning  wave  motion. 
Therefore,  the  first  tangential  mode  alone  may  exhibit  a  spinning  mode  of 
oscillations  if  the  linear  coupling  due  to  the  combustion  response  is  considered. 

Next,  check  the  condition  for  which  the  limit  cycles  exist.  From  Eqs.  (4.27) 
and  (4.28) 

ai=(-l)n*i4  (n  =  1,2,3*  •  •)  (4.34) 

This  is  the  necessary  and  sufficient  condition  for  the  existence  of  limit  cycles.  The 
remaining  task  is  to  determine  the  conditions  for  which  the  limit  cycle  is  stable. 
The  procedure  consists  of  examining  small  perturbations  in  the  vicinity  of  the  limit 
cycle.  Substituting  the  relation  rn  =  rn0  +  r'n,  X  =  Xq  +  X' ,  into  the  Eqs.  (4.20)- 
(4.22)  and  neglecting  the  higher  order  terms  yield  a  system  of  linearized  equations. 
The  solutions  are  assumed  to  vary  exponentially  in  time  and  may  be  written  in  the 
form,  r'j  =  fiexp(At),  etc.  This  leads  to  a  system  of  linear  algebraic  equations  for 
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r„,  and  X  in  matrix  form. 

^ -A  (-l)n+1014  (-l)n+1ai4r10 

(-l)n+10i4  aj  -A  (-l)nax4rio 


(-1) 


1  ^n  +  l  2on 
no 


’  n  ' 

f"4 

.  X  , 

=  0 


(4.35) 


For  the  non-trivial  solutions,  the  determinant  of  the  3x3  square  matrix  in  Eq.  (4.35) 
is  equal  to  zero,  giving  the  characteristic  equation  for  A. 


P( A)  =  A3  -  4a,  A2  +  [50?,  +  (-1)"+10?,  +  a?,]  A  =  0 


(4.36) 


Application  of  Routh-Hurwitz  criteria  gives  the  necessary  and  sufficient  condition 
for  the  stable  limit  cycle. 

ax  <  0  (4.37) 

Figure  4.1  shows  the  limit  cycle  of  the  first  tangential  mode  with  the  condition 
(4.37)  satisfied. 

Let’s  examine  this  case  from  the  acoustic  energy  conservation  point  of  view  at 
the  limit  cycle.  Acoustic  energy  density  is  defined  as 

(4'38) 

The  total  time-averaged  energy  of  the  combustion  chamber  is  obtained  using 
Eqs.  (4.1)  and  (4.13). 

£„  =  constant  •  r2n  (4.39) 


Then,  the  rate  of  energy  change  of  the  system  is 


d£n 

dt 


=  constant  •  rn 


drn 

dt 


(4.40) 


NON-DIMENSIONAL  TIME 


Figure  4.1  An  Example  Showing  the  Existence  of  Stable 


Limit  Cycle,  IT  Mode 


83 


Substituting  Eqs.  (4.20)  and  (4.21)  into  Eq.  (4.40)  and  adding  two  equations  using 
the  relation  (4.24)  yield 

^  [Si  +  S4\  =  oq(ri  +  r\)  -  2 014rir4  sinX  (4.41) 

In  the  limit  cycle,  X  approaches  Xo,  therefore,  the  right  hand  side  of  Eq.  (4.41) 
vanishes  due  to  the  condition  (4.34),  which  means  that  the  acoustic  energy  is 
conserved.  Since  aq  is  negative,  the  first  term  of  Eq.  (4.41)  represents  the  energy 
dissipated.  And  then,  energy  is  supplied  by  the  second  term,  which  represents  the 
linear  coupling  mechanism.  Careful  examination  of  linear  coupling  coefficient  014 
reveals  that  the  second  term  in  Eq.  (4.41)  is  closely  related  to  the  tangential  velocity 
effect  by  the  Eqs.  (4.106),  (4.14)  and  (4.19).  This  means  that  the  velocity- coupled 
response  function  (especially  in  the  tangential  direction)  plays  an  essential  role  for 
the  spinning  mode  of  oscillations.  This  observation  is  consistent  with  Reardon’s 
experiment35  whereby  the  velocity  effect  has  the  strong  destabilizing  effects  on  the 
spinning  mode  of  transverse  oscillations. 

4.2.2  First  Tangential/First  Radial  Modes 

Next,  the  combination  of  the  first  tangential  and  first  radial  modes  is  studied. 
The  mode  shapes  rpi,  t/>4  and  xp2  are  taken  into  account.  The  equations  governing 
the  wave  amplitudes  and  phases  of  each  mode  can  be  given  by, 

First  Tangential  Mode 

dr\ 
dt 


airj  -f  airir2  cosX 
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+a14r4  cos (fa  -fa)-  9\ 4r4  sin(<£i  -  fa) 


(4.42) 


dr  4 
dt 


=  o4r4  +  air2r4  cos  Y 

4-a^irj  cos(fa  -  fa)  +  041ri  sin(<£i  -  fa) 


(4.43) 


First  Radial  Mode 

dr2 

dt 


—  a2r2  +  bi  [rj  cos  X  +  r\  cos  Y ] 


(4.44) 


and 


dX 

dt 


dY 

dt 


20i  —  02  +  ^i2  —  2 a\r2  sinX - -  f r?  sinX  +  r?  sinKl 

r2  L  J 

2r4 

~  [#i4  cos(^i  -  fa)  +  a14  sin(^i  -  fa)] 

2 04  — 02  +  fli2  —  2<*i  r2  sin  Y  —  —  [r  ?  sin  X  +  r?  sin  yl 

r2  L  J 

2rj 


r4 


[04 1  cos(fa  -  fa)  -  a4j  sin(<^i  -  fa)] 


(4.45) 


(4.46) 


where, 


X(t)  —  2<f>i  —  fa  +  fi42^ 


Y(t)  —  2fa  —  fa  +  l)i2^ 

n12  ~  2 Wj  —  u>2 


(4.47a) 


-  fa  =  (X  -  r)/2 


And  the  same  values  of  a\ ,  a2,  6i  and  62  as  those  of  Ref.  47  are  used  in  this  chapter 


as. 


ai  =  0.1570  (|),  a2  =  -0.0521  (|) 
h  =  -0.1054(|),  62=  0.1873(|) 


(4.47  b) 
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The  following  relations  are  valid  for  the  limit  cycle. 


v\  =  Va 


XQ  —  2£i  —  £2 


(4.48) 


Yo  =  2U  -  (a 

The  reasonable  assumption  (4.30)  is  used  again;  hence,  at  the  limit  cycle, 
Eqs.  (4.42)-(4.46)  become 

0  =  Ox  +  air2o  cos X0  +  ox4  cos  ^  -  614  sin  (4.49) 

„  .  „  ,  (*o-lo)  ,  „  .  (X0-Y0) 

0  =  a4  +  air2o  cos  Yo  +  <*4i  cos - - - +  041  sin - - - -  (4.50) 


0  =  o2r20  +  feir^o  [cos  XQ  +  cos  y0] 


(4.51) 


0  =  —26i  — 02 +  fii2  —  2aiT2o  sinXo - — r?0  [sin X0  +  sin  Y0] 

r20 


2 r40  „  (X0  -Y0)  ,  .  (X0  -Yo) 

-  OlA  COS - - - +  oi4  sin - - - 

rao  Z  Z 


0  =  -204  -02  +  fii2  -  2axr20  sin  Y'o - -r?0  [sin.Y0  +  sin  Y'o] 

r2  1  1 


(4.52) 


2r,„  „  (Xn-Y„)  .  (X0-Y„) 

-  flu  cos - - - 041  sm - s - 

t4o  *  Z 


(4.53) 


Spinning  wave  motion  exists  only  when  the  phase  difference  should  be  t2n+1)"\ 
FVom  Eqs.  (4.49)-(4.53),  only  the  condition  r2o  — »  0  makes  the  system  a  spinning 
mode.  When  r2o  approaches  0, 


{,=«<  + ^4^  (n  —  1,2,3*..) 


(4.54) 
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Now  check  the  condition  for  which  the  limit  cycles  exist.  From  Eqs.  (4.49)  and 
(4.50),  the  same  condition  as  (4.34)  in  Section  4.2.1  is  obtained. 

ai=(-l)n014  (n  =  1,2,3-  •  •)  (4.55) 

This  is  the  necessary  and  sufficient  condition  for  the  existence  of  limit  cycles.  The 
remaining  task  is  to  determine  the  conditions  for  which  the  limit  cycle  is  stable. 
However,  the  peculiar  structure  of  the  r2o  equation  does  not  make  it  possible  to  get 
the  conditions  for  which  the  limit  cycle  is  stable.  Figure  4.2  shows  the  limit  cycle  of 
the  1T/1R  modes,  in  which  the  pattern  r2o  — *■  0  can  obviously  be  seen  at  the  limit 
cycle. 

Next,  check  this  case  from  the  energy  conservation  point  of  view.  From 
Eqs.  (4.42)-(4.44),  the  rate  of  energy  change  for  the  1T/1R  modes  becomes 
d 

—  [Si  +S4  +  £2]  =  «i(rj  +  rl)  -  2^14r1r4  sin -  <j>4) 

(4.56) 

+  a2rj  +  bir2(rl  cosX  +  r\  cosF) 

It  can  easily  be  shown  that  the  right  hand  side  of  Eq.  (4.56)  vanishes  at  the  limit 
cycle  by  substituting  condition  (4.55).  Therefore,  the  conclusion  is  reached  that 
there  is  no  contribution  of  the  first  radial  mode  to  the  spinning  mode  of  oscillations 
after  the  limit  cycle  is  attained  for  the  combination  of  1T/1R  modes.  That  is,  the 
contribution  of  the  first  radial  mode  is  limited  in  the  initial  stages,  but  its  effect 
disappears  once  the  limit  cycle  is  obtained,  as  can  be  clearly  seen  in  Fig.  4.2. 
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4.2.3  First  Tangential/Second  Tangential  Modes 

Mode  shapes  ip j,  ip 4,  ipz  and  ip 5  are  considered  for  the  combination  of  the 
first  and  second  tangential  modes.  The  analysis  of  1T/2T  modes  follows  the  same 
approach  as  that  of  the  1T/1R  modes  for  the  spinning  mode  of  oscillations.  The 
following  relations  for  the  linear  parameters  and  linear  coupling  parameters  are 
already  found  in  Section  4.1. 


<*1  =  <*4,  #1  =  04 

<*3  =  q5,  03  =  05  (4.57a) 

<*14  =  <*41 1  014  =  041 

<*35  =  ~ <*53>  035  =  “053  (4.576) 


The  approximate  equations  governing  the  amplitudes  and  phases  of  each  mode  are 
given  by, 


First  Tangential  Mode 


dr\ 

dt 


<*iri  +  a2rir3  cosX  +  a2r4r5  cos  Z 


+  a14r4  cos(</>i  -  <p4)  ~  014^4  sin(<^  -  <p4) 


dr4 

~dt 


a4r4  —  a2r3r4  cosF  +  02^5  cosZ 


+  a4iri  cos(^i  -  <p4)  +  04irj  sin(<£i  -  <p4) 


(4.58) 


(4.59) 


Second  Tangential  Mode 
dr3 


dt 


=  <*3r3  +  62  [rj  cos  X  -  r4  cos  F] 
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+  a35r5  cos(^3  -  fa)  -  035r5  sin(<£3  -  fa) 

dr  5 

—  =  a5r5  +  2b2rir4  cos Z 
at 

+  a53r3  cos(fa  -  fa)  +  @53r3  sin(fa  -  fa) 


and 


dX 

dt 


—  —  2#i  +  $3  +  fii3 - 

r3 


2a2r4r5  . 


■  „  2  r4 

sinZ - 

ri 


[2a2rf  +  b2r\\  sin  X  +  siny 

r3 

[<*14  sin(<^>i  -  fa)  +  014  cos(^!  -  fa)] 


7*5 

+  —  [<*35  Sin(^3  -  <£5)  +  #35  COs(fa  -  fa)] 
r  3 


dY 

dt 


- -  5-  sin  Z  +  [a4j  sin(<£i  —  fa)  —  041  cos(<fn  —  <j)4)\ 

T5 

H - [<*53  sin(<^3  -  fa)  +  053  cos(<£3  -  <£5)] 

T3 


dZ 

di 


=  —  &i  —  04  +  65  +  fi13  —  a2r3  [sinX  —  sinF] 


-  02 


/r4r5  rirs  \  262r!r4 

v  ri  ta  / 


r  5 


sin  Z 


V\ 

—  [a44  sin(^i  -  <£4)  +  0i4  cos(<^i  -  fa)] 
r  i 


7*2 

+  —  [<*4i  sin(fa  -  fa)  -  641  cos(fa  -  fa)] 
r4 


- [<*53  sin(^3  -  fa)  -  053  cos(fa  -  fa)] 

r5 


(4.60) 

(4.61) 

(4.62) 

(4.63) 


(4.64) 
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where, 

X  =  2(^i  —  (f)  3  +  Q\3 1 
Y  =  2<f>4  —  <f>3  +  f2i3< 

Z  —  <f>\  +  <f>4  —  <f>5  +  fll3^ 

(4.65) 

fii3  =  2a.’ j  —  u>3 

fa  -  fa  =  (X  -  Y)/ 2 
<6  -  fa  =  Z  -  (X  +  Y)/2 

At  the  limit  cycle,  the  same  procedure  as  that  in  Section  4.2.2  is  followed  to  find 

i'l  =  ^4,  V3  =  i/5 

*o  =  26  —  6 

(4.66) 

Yq  —  2^  —  £3 
Zq  —  —  6 

The  Following  assumption  for  the  limit  cycle  is  used,  which  is  practical  in  a 
reasonable  sense. 

rio  =  r40,  r30  =  r50  (4.67) 

Considerable  manipulations  of  Eqs.  (4.58)-(4.64)  yield  the  following  relationship. 


X0  =  Y0  +  (2n  +  1)tt  (n  =  1,2,3-  •  •) 
( Xo+Yo )  ,  (2*  +  1)tT 


Z0 


+ 


2  2 

Then,  the  relations  between  frequency  modulations  are 

(2n  +  1)tt 


(*  =  1,2,3---) 


6=6  + 

6  =  6  + 


2 

(2*  -f  1  )tt 


(n  =  1,2,3---) 
(Ar  =  1,2,3  -  •) 


(4.68a) 

(4.686) 

(4.69a) 
(h  J96) 
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This  means  that  the  phase  differences  between  mode  1  and  mode  4,  and  between 
mode  3  and  mode  5  are  f  or  ^  at  the  limit  cycle.  On  t^e  other  hand,  acoustic 
pressure  for  the  combination  of  the  first  and  second  tangential  modes  is 
p'(r,t)/p  =  r10  sin^f  4-  <^iO)cos0J1(/c1r) 

4-  r40  sin(uqt  +  ^40) sinflJ^/cjr) 


+  r30  sin(u>3t  +  <f>30)  cos  20J2(/c3r) 


(4.70) 


4-  r50  sin(a,'3t  +  ^50)  sin 2^  J2(K3r) 

Eq.  (4.70)  can  be  further  simplified  using  Eqs.  (4.67)  and  (4.69)  to  give 
p\r,t)/p  =  rjo  sin[(u>!  4-  Ui )<  +  6  -f-  «7i(/c1r-) 

(4.71) 

+  r3o  sin [(u>3  4-  v3 )t  4-  £3  4-  20]  J2(k 3r ) 

Equation  (4.71)  shows  that  the  spatial-  and  time-dependent  paxts  of  the  first  and 
second  tangential  modes  are  combined  together,  resulting  in  the  spinning  wave 
motion.  Therefore,  the  combination  of  1T/2T  modes  exhibits  the  spinning  mode 
of  oscillations  within  the  chamber  if  the  linear  coupling,  due  to  the  combustion 
response,  is  taken  into  account.  Some  manipulations  of  the  Eqs.  (4.58)-(4.61)  give 
the  solution  for  the  amplitude  of  each  mode  at  the  limit  cycle. 


rio  — 


T30 


[Ql  —  (~ l)n014]  [c*3  —  (  — 1)^35] 

4a262  cos2  Xo 
Ql  ~  (-l)n0!4 


(4.72) 

(4.73) 


2a2  cos  Xo 

According  to  the  definitions  of  a2  and  b2  in  Eq.  (4.47),  a2b2  is  negative;  hence,  for 


the  real  value  of  r^,  the  following  relation  should  be  satisfied. 

[a,  -  (-!)"«»]  [a,  -  (-1)*93»]  <0  (n,k  =  1,2,3  •  ■  •) 


(4.74) 
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This  is  the  necessary  and  sufficient  condition  for  the  existence  of  the  limit  cycle. 

The  next  step  is  to  determine  the  conditions  for  which  the  limit  cycle  is  stable. 
Substituting  the  relations  rn  =  r„o  +  r'n,  X  =  Xo  +  X' ,  etc.  into  the  Eqs.  (4.58)- 
(4.64)  and  neglecting  higher-order  terms,  the  system  of  the  linearized  equations  is 
given  by. 

dr[  i  [<*1  +  (— 1)”014  ,  ,  |"—  <*1  +  (  —  l)n#i4 

*■  =  r>[ — 2 — J+r<[ - 2 - j 

+r'3  (a2r10  cosX0)  +  r'5  ( a2r10  cosX0) 

+X '  -a2raor3osinXo  +  (-l)n+1^y^ 

+r  (_l)„£^io  +Z'(-a2r10r3oSinXo)  (4.75) 

dr4  ,  [—  ai  +  (  —  1  )”^i4 1  ,  _*  [oi  +  (  — l)n0i4 

*■=  r>[  2  j+rn  2  j 

+r3  (02^10  cosXo)  +  r'5  (a2r10  cosX0) 

+X'  (  — l)n — - —  +  Z'  (— a2rjor3o  sin  Xo) 

+Y'  -a2r10r30sinX0+(-l)n+1^yi2-  (4.76) 

~  =  r\  (2fe2rio  cosXo)  4-  r4  (2 &2r10  cosX0)  +  r3(a3)r5  ((-1)*+1035) 

+X'  -b2r\Q  sinX0 +  (-!)* 

+r  -&2r?0sinX0  +  (-l)*^|^  +  Z'[(-l)"+1a35r30]  (4.77) 

~  =  rj  (262rxo  cosXo)  +  r4  (262rio  cosXo)  +  r'3  ((  — l)fc+1^35)  +  r'5(a3) 
at 
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+A' 


(-1) 


k+ 1  <*357\jO 


+  Y' 


+  ('l) 


k+1  ^35^30 


-*-Z'  [— 2fc2rio sinAo  +  (  —  l)*035no] 


and 


dX' 

dt 


=  ri 


/  2fl2^30  .  v"  262^10  v1  ,  (  1  \n 

,  - smXo - smXo  +  (-l) 

no  no 

— 2a2no  .  v 
- sin  Aq  — 


no 

'  ■  Y  2 

—2a2  sin  Ao  H - ; 

r; 


^i-0  sin  A0 -(-!)' 


30 


:  Q35 
30  J 


+K 

+^3 

+^5 

+X' 

+Y' 

+  Z'  [— 2a2r30  COS  Ao  -  (-1)*^35] 


2&2n0  .  v  /  ,n2o14 

- sin  A0  —  (  — 1)  — 

no  r 


2qm 

no  J 

2'>i, 
r10 


ka  35 

no  J 


— 2a2  sin  Ao  +  +(  — l)1 

-2a2r3o  cos  Ac  -  hlkcosx0  +  (-l)"014  +(-l)* 
no 

Z^IO  cosXo  _  (_1}n^4  +  (-1)*^ 

no  ^ 


dr 

dt  ~  ri 


+n 


— 2a2no  .  v  262^10  v-  ,  ,  1V 

- sin  A0 - sin  A0  +  (  —  1) 

no  no 

2a2n0  ■  v  262^10  »  v  1  f  1  \n 

- sinAo - sin  Ao  +  ' -  lj 

no  no 

,2 


+7-3 

+^5 

+A' 


1V>2°14  1 

no 

J 

,n  2ai4 

1  r,„J 

ka  35 


— 2a2  sinAo  +  —a  -1-0  sinAo  —  (  —  1) 

*30  r30 


— 2a2  sinAo  d-  d"( —  1)* 


035 

30 


cos  Ao  -  (  1  )n $14  +  (-1)*% 

no  z 
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+Y' 


—2a2r30  cos  Ao  — 


V|o 

r30 


cosJf0  +  (-l)”«i4  +(-!)*% 


* \-Z '  [—  2a2T3Q  COS  Xq  —  (  —  l)*$3s] 


(4. SO) 


dZ' 

dt  ~  ri 
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+X' 

+Y' 
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sin  Ao 
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2hrio  .  v 
sin  Ao 


r30 


— 2g2  sin  Ao  "h  +( —  1)* - 


30 


—2a2  sin  Ao  -f 


2b2r^0  . 


*  V  /  -1  \k  ^35 

smA0  -  (-1)  - 


30 


r30 


ir  /  ,  ,fc^35 

-a2r30  cos  A0  -  (-1)  — 

ir  /  i\fc^35 

-a2r30  cos  A0  -  (-1)  — 


-2a2r30  cosX0  —  -~?riQ-  cos  A0  +  (-l)fc+1035 


**30 


(4.81) 


If  these  equations  are  expressed  in  matrix  form,  it  gives  a  7x7  matrix  which  is 
impossible  to  handle  analytically.  For  the  analytical  tractability,  let  us  take  the 
assumption  that  r[  —  r'A  and  r'3  —  r'5,  and  X '  —  Y'  —  Z'  which  means  the  initial 
phase  differences  between  mode  1  and  mode  4,  and  between  mode  3  and  mode  5 


are  |  or  y.  The  system  is  then  reduced  to  a  3x3  matrix  that  can  be  handled 


analytically;  then,  the  characteristic  equation  for  A  is  given  by, 


where, 


A3  -f  Pi  A2  P2\  P3  —  0 

Pl  =  _2 (ai  +  a' ) 

P2  =  a3  [03  -  (4a',  -  03)  tan2  A0] 

P3  =  2a',  a!,  [(2a',  +  a'  )(1  +  tan2  A0)] 


(4.82) 


(4.83) 
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Application  of  Routh-Hurwitz  stability  analysis  yields  the  following  conditions  for 
the  stable  limit  cycles. 


Ot\  4-  c*3  <  0 


2a  j  +  03  <  0 

,  Hk|  /-X2  +  2x  +~2 
1  2  +  xy  X2  +  2x  +  2 


if  <  0 


(4.84) 


where, 

Ot\  —  Qi  -  (-1)"^14 
03  =  a3  -  (-l)fc^35 


k  =  2  (61  +  (— l)"ai4)  —  (63  +  (  —  1)^035)  —  (2u>\  —  u>3) 


(4.85) 


The  linear  coupling  parameters  are  incorporated  with  the  linear  parameters.  Figure 
4.3  shows  the  overall  result  of  Eq.  (4.84)  for  the  special  case  of  k=0.025,  where  the 
shaded  regions  represent  the  conditions  for  the  stable  limit  cycles. 

Figure  4.4  shows  the  time  history  of  the  limiting  amplitudes  of  1T/2T  modes  with 
parameters  which  satisfy  the  conditions  in  Fig.  4.3.  The  stable  limit  cycles  can  be 
clearly  seen. 

4.2.4  First  Tangential/First  Radial/Second  Tangential  Modes 


Spinning  oscillations  of  the  combination  of  the  first  three  modes  are  also 
investigated  using  the  same  analysis  as  that  of  the  two-mode  combination  cases. 
The  approximate  equations  governing  the  amplitudes  and  phases  of  each  mode  are 
as  follows, 


NON-DIMENSIONAL  TIME 


Figure  4.4  An  Example  Showing  the  Existence  of  Stable 


Limit  Cycle,  1T/2T  Modes 
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First  Tangential  Mode 


=  a\r i  +  airir2  cosX  +  a2r\rz  cosF  +  a2r4r5  cos Z 


+  ai4r4  cos(^i  -  (f>4)  -  0i4r4  sin(^»i  -  <f>4) 


(4.86) 


=  a4r4  +  a4r2r4  cos X\  —  a2r3r4  cos  Y\  +  a2r\r3  cos Z 
at 

+  or4in  cosfOj  -  04)  +  041rj  sin(0)  -  04  j 


(4.87) 


a<p  i 

—r~  —  —  B\r\  —  a\r\r2  sinX  —  a2r4r3  sinF  —  a2r4r5  sin  Z 
at 

-  ai4r4  sin(^>i  -  <j>4)  -  044r4  cos (<f>x  -  <f>4) 


(4.88) 


\-jj-  =  —  04r4  —  a4r2r4  sinXi  +  a2r3r4  sinFj  —  a2rir5  sin Z 
+  a4iri  sin(^i  -  <f>4 )  -  041ri  cos^j  -  <f>4) 


(4.89) 


First  Radial  Mode 


=  a2r2  +  b\  [rj  cosX  +  r\  cosXj] 


(4.90) 


!  —j7~  =  ~  B2r2  +  bi  [r\  sinX  +  r\  sinXi] 


(4.91) 


Second  Tangential  Mode 


=  a3r3  -f  b2  [rj  cos  F  -  r4  cos  Fi] 

+  a35r5  cos(^3  -  <f>b)  -  035r5  sin(<£3  -  fc) 


(4.92) 


~rr  =  a5r5  +  2b2ri  r4  cos  Z 
at 

+  0:53^3  cos(<^3  -  <t>5)  +  033r3  sin(^3  -  4>s ) 


(4.93) 
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r3  —  =  -  63r3  +  b2  [rj  sin  Y  -  r\  sin  Yi] 

-  «35^5  sin(^»3  -  <f>5)  -  035 r5  cos(<^3  -  <f>5 )  (4.94) 

r5  =  -  05r5  +  2b2rir4  sin  Z 
at 

+  ^53^3  cos(<^>3  -  <t> 5)  -  053r3  sin(^>3  -  <f>5 )  (4.95) 


where, 

X  —  2(f>\  —  <f>2  +  £li2t 
Y  =  2<f>i  —  <f> 3  4-  fij3 1 
X\  —  2(f) 4  —  <f>  2  +  Oi2< 
Yi  =  2<f>4  —  <f)  3  +  ftl3^ 


(4.96) 


Z  —  <f>\  +  <f>4  —  (f>$  +  J^i3  / 

There  are  ten  independent  equations  for  ten  unknowns,  but  the  algebra  is  much 
more  complicated  than  that  for  the  two-mode  combination  cases.  Direct  numerical 
integration  of  Eqs.  (4.86)-(4.95)  is  performed  with  various  combinations  of  linear 
and  linear  coupling  parameters.  Similar  to  the  result  of  the  1T/1R  modes,  the 
amplitude  of  the  first  radial  mode  r2o  goes  to  zero  for  the  spinning  wave  motions. 
Figure  4.5  shows  the  time  history  of  the  amplitudes  for  the  1T/1R/2T  modes  case, 
in  which  the  pattern  r2o  — >  0  can  be  clearly  seen.  Therefore,  the  first  radial  mode 
in  the  combination  of  the  1T/1R/2T  modes  has  no  contribution  after  the  limit, 
cycle  is  reached.  But  in  the  initial  stages,  it  helps  the  system  reach  the  limit  cycle 
regardless  of  the  initial  phase  differences  and  once  the  limit  cycle  is  attained,  its 
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effect  disappears.  The  phase  difference  of  the  first  tangential  modes  and  that  of 
the  second  tangential  modes  always  become  odd  number  multiples  of  That  is, 
the  system  always  degenerates  to  the  spinning  wave  motion  regardless  of  the  initial 
phase  differences. 


PRESSURE  AMPLITUDE 
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CHAPTER  5 

NONLINEAR  ACOUSTIC  WAVES  IN  BAFFLED  CHAMBERS 

This  task  represents  an  outgrowth  of  the  linear  acoustic  analysis  described 
in  Chapter  3.  A  linear  acoustic  analysis  provides  a  means  for  assessing  whether 
a  particular  system  is  stable  or  unstable  with  respect  to  small  perturbations,  but 
gives  no  information  about  the  ultimate  amplitudes  that  unstable  motions  reach.  In 
practical  systems,  constant  amplitude  oscillations  (i.e.  limit  cycles)  have  often  been 
observed.  Since  the  occurrence  of  a  limit  cycle  in  a  dynamical  system  requires  the 
influence  of  nonlinear  interactions,  modeling  of  this  phenomenon  within  a  chamber 
requires  appropriate  consideration  of  the  nonlinear  terms  in  the  governing  equations. 
There  have  been  many  studies47 ,48,56  addressing  nonlinear  analysis  in  the  perfect 
(unbaffled)  cylindrical  chambers  but  none  exists  for  the  baffled  combustors.  This 
chapter  deals  with  the  nonlinear  acoustic  analysis  in  baffled  environments. 

5.1  Nonlinear  Acoustic  Oscillations 

The  wave  equation  derived  in  Chapter  2  contains  various  nonlinear  processes 
of  concern.  Among  all  nonlinear  processes,  the  effect  of  second  order  nonlinear  gas 
dynamics  is  examined  systematically  here. 

Unlike  the  linear  acoustic  analysis,  the  simple  assumption  of  time-harmonic 
behavior  of  all  the  acoustic  variables  is  no  longer  valid.  However,  for  a  typical  rocket 
engine,  the  mean-flow  Mach  number  and  the  variation  of  mean  pressure  are  small. 
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Under  such  conditions,  the  frequencies  and  the  spatial  variation  of  the  unsteady 
motions  deviate  only  slightly  from  the  classical  acoustic  field  obtained  for  the  same 
geometry,  but  without  any  source  terms  (h  =  n  •  f=0).  Thus,  the  solution  to  the 
wave  equation  can  be  faithfully  approximated  by  a  synthesis  of  the  linear  normal 
acoustic  modes  pn( r),  with  unknown  time- varying  amplitudes  r)n(t )  as  coefficients: 

P'=P$>  «(*)?»(*)  (5-1) 

n 

Correspondingly,  the  acoustic  velocity  field  is  written  as 

«M)  =  (5-2) 

where  the  linear  mode  shape  pn(r)  satisfies  the  Helmholtz  equation  throughout  the 
chamber,  subject  to  rigid  surface  boundary  conditions. 

V2pn  +  klpn  —  0  (5.3a) 

n  •  Vpn  =  0  (5.36) 

The  time- varying  amplitudes  rjn(t )  give  the  temporal  nonlinear  behavior  of  the 
system.  Eqs.  (5.1)  and  (5.2)  are  not  exact  representations  of  the  true  fields,  for  the 
boundary  conditions  are  not  satisfied;  Eq.  (5.1)  gives  n  •  Vp'  =  0  because  the  p„ 
satisfies  Eq.  (5.36)  and  according  to  Eq.  (5.2),  the  velocity  fluctuation  vanishes  on 
the  boundary.  This  means  that  Eqs.  (5.1)  and  (5.2)  do  not  accurately  reproduce 
the  spatial  structure  of  the  unsteady  motions  near  the  boundary.  Nevertheless,  the 
errors  are  small  when  the  perturbations  contained  in  h  and  f  of  Eqs.  (2.30)  and 
(2.32)  are  small.  Furthermore,  due  to  subsequent  spatial  averaging,  the  equations 
found  for  the  amplitudes  will  provide  a  satisfactory  basis  for  studying  real  problems. 
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Note  that,  owing  to  the  geometric  discontinuity  set  by  the  baffles,  the  linear 
mode  shape  pn  is  not  the  same  as  its  counterpart  for  a  perfect  cylinder.  Rather,  it 
is  expressed  as  a  series  of  transverse  normal  modes. 

Pn(r)  =  £Cni(z)Mrf0)  (5.4) 

i 

Here,  ipi(r,0)  is  the  transverse  normal  mode  function  for  a  cylinder.  The  axial 
distribution  function  £„,(x)  can  be  determined  using  the  eigenfunction-expansion 
technique  discussed  in  Chapter  3.  Thus,  the  nonlinear  solution  found  in  this  chapter 
is  constructed  using  the  linear  acoustic  solutions  found  in  Chapter  3. 

A  system  of  ordinary  differential  equations  in  terms  of  the  time-varying 
amplitudes  of  Eq.  (5.1)  can  be  obtained.  This  is  done  by  averaging  the  conservation 
equations  over  the  volume  of  the  chamber  using  a  version  of  the  method  of  least 
residuals,  which  is  equivalent  to  Galerkin’s  method  with  the  linear  mode  shape  as 
a  weighting  function.  Now  multiply  Eq.  (2.29)  by  pn  and  Eq.  (5.3a)  by  p',  subtract 
the  results,  integrate  over  the  entire  volume,  and  apply  Green’s  theorem  to  get 

+  k"p"p')dv  =~JII  p"hdv  ~  §P"fnds  (»•*) 

Substitution  of  Eqs.  (5.1)  and  (5.4)  into  Eq.  (5.5)  and  some  rearrangement  of 
the  result  leads  to  a  system  of  second-order  ordinary  differential  equations  for  the 
temporal  behavior  of  nonlinear  oscillations. 

^  4"  ~  (5-6) 

m 

where  the  amplitude  coefficients,  Cnm ,  and  forcing  function,  Fn,  are  defined  as, 

Cnm  =  Y,E>  j 


(5.7) 
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pnhdV  +  j^>pnf  ‘ 

where 

Ei=  JJ  tf(r,8)rdrd0  (5.9) 

Equation  (5.6)  represents  the  time-dependent  motions  of  a  collection  of  nonlinear 
oscillators,  one  oscillator  being  associated  with  each  linear  mode.  The  principal 
difference  between  the  present  analysis  and  the  unbaffled  cases47,48  lies  in  the 
structure  of  the  ordinary  differential  equation  (5.6),  in  which  strong  coupling  among 
modes  occurs  on  both  sides  of  the  equation.  Now,  the  problem  is  reduced  to  solving 
this  system  of  second-order  ordinary  differential  equations  for  the  time-varying 
amplitude,  given  a  specified  initial  condition. 

The  forcing  function,  Fn  represents  the  deviation  from  unperturbed  acoustic 
modes  and,  therefore,  provides  the  mechanism  that  allows  for  initial  finite 
oscillations  to  grow  or  decay.  As  a  result,  the  various  components  of  Fn  are  essential 
to  the  study  of  combustion  instabilities  in  the  baffled  combustion  chambers.  It  is 
convenient  to  classify  these  components  of  Fn  according  to  the  orders  of  p'  and 
in  each  term.  Substitution  of  Eqs.  (2.30)  and  (2.32)  into  Eq.  (5.8)  yields  the  explicit 
expression  of  the  forcing  function  Fn.  After  some  straightforward  manipulations, 
Fn  can  be  written  as 


M  [d\ 
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l«]  W 


+  (^Pn^(Ugp'  +  U gP  +  Ug'p ')  *  II  <fS 

S -  '  - - - 

W 

-///(^•vP»)  +  (ip„  fV 

V|  ■  ■——■  ■■  ■  ^  y  - 

(i) 

The  first  four  terms,  [a]  —  [d],  axe  referred  to  as  the  linear  acoustics,  and  the 
next  four  terms,  [e]  —  [h\,  the  second  order  nonlinear  acoustics.  The  last  two  terms, 
[t]  —  [7],  are  related  to  the  combustion  response  and  other  sources  such  as  particle 
damping.  If  only  the  nonlinear  contribution  from  the  gasdynamics  of  the  unsteady 
flow  is  taken  into  account  in  this  study,  the  forcing  function  can  be  expressed  in 
terms  of  the  wave  amplitude  r)n,  with  the  form 

OO  OO  OO 

En  =  ^  ,  [Dni’ni  d-  EniT) j]  ^  ^  ^  ^  \_^nijViVj  d*  Bnij'r)i'rlj\  (5.11) 

i=l  t=l  >= 1 

Note  that  the  cross  coupling  term  r\r\  does  not  appear  in  the  formulation 
of  Eq.  (5.11)  since  nonlinear  combustion  responses  and  the  interactions  between 
mean  flow  and  nonlinear  acoustics  are  not  included  in  this  analysis.  The  explicit 
forms  of  the  above  coefficients  are  given  by, 

Dni  =  IIIl11  +  “  T*(vpi  x  (v  *  a,)) ' Vi‘" 
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+  IfPnPi  V  •  U g 


dV 


(5.12a) 


and 


(5.126) 


Afiij  — 


Trrrri  (**2  +  */)(2*.2  +  2 k?  -  kn2) 

S'rki  kj  L 

4  oz.L.2l  2 


'n«j 


_  a2 Ini] 


dnij  — 


47 

JJJ  PnPiPj  dV. 


—  kn  ~  Byki  kj 


2(ki2  +  kj2)(l-i)-kn' 


(5.12c) 

(5.12d) 

(5.12c) 


The  linear  coefficients  Dn,  and  Eni  arise  from  all  linear  processes  in  Eq.  (5.10), 
and  the  nonlinear  coefficients  Anij  and  Bnij  represent  the  second  order  nonlinear 
acoustics  part  of  Eq.  (5.10).  Note  that  all  coefficients  depend  on  the  linear  mode 
shapes  and  frequencies  of  the  chamber. 

After  the  forcing  function  F„,  Eq.  (5.11)  is  substituted  in  Eq.  (5.6),  this  can 
be  rewritten  in  a  matrix  form: 


where 


[c]w  =  m 


[C]{rj}=  Y,0™’)'" 

m=l 

oo 

{/}  =  -  +  (En,  +  Crivl)Vi 

1=1 


oo  oo 

-EE*  nij^li^lj  "1”  ^nijViVj 

.=  1  j=  1 


(5.13) 


(5.14a) 


(5  1  tb) 
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Now  the  problem  is  reduced  to  solving  a  set  of  simultaneous  differential  equations 
(5.13).  But  unlike  the  unbaffled  case,  the  wave  motion  can  no  longer  be  viewed 
as  a  weakly  coupled  system.  Consequently,  the  method  of  time-averaging  used  by 
Culick,46  which  can  replace  a  system  of  second-order  ordinary  differential  equations 
by  an  equivalent  system  of  first-order  equations,  cannot  be  practically  applied  to 
solving  Eq.  (5.13)  due  to  this  strong  coupling  among  modes.  Therefore,  numerical 
integration  of  Eq.  (5.13)  is  performed  using  a  fourth-order  Runge-Kutta  scheme  to 
get  the  time  history  of  the  amplitudes  r/n(t).  However,  a  similar  format  as  used  in 
the  method  of  time  averaging  is  employed  for  the  representation  of  the  amplitude 
of  the  limit  cycle.  Namely,  the  time-varying  amplitudes,  r/n(<)  are  transformed  to 
two  variables,  rn(t)  and  <f>n(t)  as  follows, 


1r.(<)  =  r„(()sin[u>„f  +  ipjt)\ 


(5.15) 


where  rn(t)  and  <f>n(t)  are  the  amplitude  and  phase  of  the  limit  cycle,  respectively. 
Then  the  amplitude  of  the  limit  cycle,  r„(f),  takes  the  form, 


rn(t)  =  \l  M*)]2  r 


Vn(t) 


u 


n  J 


(5.16) 


This  expression  is  used  to  represent  the  temporal  variation  of  the  amplitude  of  the 
limit  cycle.  The  major  task  at  this  point  is  to  determine  the  limiting  amplitude 
rn(t )  characterizing  the  nonlinear  oscillations. 

For  the  two-dimensional  rectangular  baffled  chamber,  the  results  from  the 
two-dimensional  linear  acoustic  analysis  are  used  to  determine  the  coefficients 
Cnm  characterizing  the  baffled  system.  In  calculating  the  coefficients  C„m,  the 
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assumption  of  no  mean  flow  is  made  to  allow  concentration  upon  the  nonlinear 
behavior  itself.  The  calculation  is  performed  in  two  regions  separately  (i.e.  baffle 
compartments  and  the  main  chamber)  and  then  combined  to  give  the  final 
expression  for  the  coefficients.  That  is, 


where 


(5.17a) 


(5.176) 


Here,  the  B  and  C  subscripts  denote  integration  in  the  baffle  compartments  and 
the  main  chamber  regions,  respectively. 

The  calculation  of  IntJ  in  Eq.  (5.12e)  is  also  treated  separately, 


lnij 


=  Jjj  PuPiPjdV 

=  ^2  5Z  f  /  CnaQbCjcdx  ■  l 

a  6  c  L./0  JE 


rl’a'Pb'Pcdy 


(5.18) 


+  /  CnaQbCjcdx  ■  /  IpaMcdy 
JLt  JC 


Equation  (5.18)  is  used  to  calculate  nonlinear  parameters  Anij  and  BniJ  from  (5.12c) 
and  (5.12 d),  respectively. 


For  the  three-dimensional  cylindrical  baffled  chamber,  the  same  procedure  as 
the  two-dimensional  chamber  case  is  followed.  Then,  the  co<  u  nt  -  Cnm  axe 


no 


determined  as 


-?[«/ 

-?K 


CniCmidl 


Cn«Cni»dx 


+ 


?K 


CniCmidx 


where 


J3?  =  JJ  il>?(r,0)rdrdO 


(5.19  a) 


(5.196) 


For  the  determination  of  nonlinear  parameters  AniJ  and  BntJ,  the  coefficient  Jnij 
in  Eq.  (5.12e)  is  calculated  first  as 


Inij 


=  JJJ  PnPiPjdV 

=  EEE[jf 

a  b  c 


A  *  A 

Cna  CibCjc^1 


IJ- 


i/>i,  V>crdrd0 


(5.20) 


+  /  c naCibCjcdx  ■  /  /  V’ftV’crdrdfl 

JLt  JJc 


5.2  Results  and  Discussion  of  Nonlinear  Acoustic  Analysis 

From  the  energy  point  of  view,  linear  oscillations  can  be  visualized  as  having 
completely  uncoupled  modes  of  energy  storage  in  the  oscillating  fields.  However,  it 
is  highly  unlikely  that  a  real  system,  certainly  a  system  as  complicated  as  a  liquid 
rocket  engine,  will  be  completely  linear  with  all  modes  uncoupled.  If  the  different 
modes  are  coupled  with  each  other,  it  is  quite  possible  that  draining  energy  from 
one  mode  will  also  drain  energy  from  other  modes  through  the  nonlinear  processes. 


Figure  5.1  Block  Diagram  Showing  the  Energy  Transfer 
Relation  among  Modes 

Therefore,  given  the  circumstances  that  instabilities  are  likely,  it  is  necessary  to  pay 
attention  to  nonlinear  behavior.  The  overall  energy  transfer  relationship  among 
modes  through  nonlinear  interactions,  and  how  this  is  related  with  linear  processes, 
are  conveniently  outlined  in  the  block  diagram  in  Fig.  5.1. 

Useful  analytical  results  from  the  nonlinear  analysis  described  in  the  previous 
section  can  fully  describe  the  time-varying  amplitudes  in  finite  oscillations  for  a 
baffled  chamber.  The  complexity  of  the  problem  depends  on  how  many  nonlinear 
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modes  fire  to  be  considered  in  the  analysis.  Unfortunately,  but  as  is  expected,  the 
analysis  becomes  more  and  more  complicated  as  the  number  of  modes  considered 
in  the  nonlinear  synthesis  is  increased.  Because  of  the  complexity  of  the  actual 
analysis,  only  the  case  of  the  first  two  modes  combination  (i.e.  n  —  1  and  n  —  2)  is 
treated  for  the  two-dimensional  analysis.  This  case  represents  the  simplest  possible 
situation  and  can  serve  as  a  basis  for  analyzing  more  complicated  problems.  Culick 
and  Yang9  showed  that  the  two  mode  expansion  is  valid  over  a  broad  range  of 
linear  and  nonlinear  parameters  that  govern  the  global  quantitative  behavior  for 
unbaffled  chambers.  Although  higher  modes  have  the  capability  to  modify  the 
stability  domain,  the  amplitudes  of  higher  modes  for  many  practiced  systems  are 
negligibly  small  due  to  the  efficient  viscous  deimping  at  high  frequencies.  Thus, 
the  description  of  nonlinear  acoustic  fields  in  the  baffled  chambers  by  the  first  two 
modes  may  be  justified. 

The  coefficients  Cnm,  which  represent  the  major  distinction  between  the  baffled 
and  the  unbaffled  chamber  nonlinear  analysis,  are  functions  of  both  baffle  length 
and  chamber  oscillation  frequency.  Although  the  combinations  using  only  the  first 
two  nonlinear  modes  are  under  examination  (n  =  1  and  n  =  2),  the  number  of 
simultaneous  equations  to  be  solved  increases  by  a  factor  of  2m,  where  m  is  the 
number  of  linear  modes  considered  in  one  nonlinear  mode.  However,  m  should  take 
on  as  many  values  as  possible.  This  makes  the  analysis  complicated  to  a  great 
extent,  since  even  a  few  modes  in  the  expansion  mean  a  large  set  of  equations  to 
be  solved.  For  this  analysis,  the  chamber  with  baffle  length  Lf,  =  0.2 H  is  used  and 
the  frequencies  obtained  from  the  linear  acoustic  analysis  are  employed  to  estimate 


113 


these  coefficients  from  Eqs.  (5.17)  and  (5.19)  for  the  two-dimensional  and  three- 
dimensioned  chambers,  respectively.  Values  of  Cnm  are  tabulated  for  the  first  three 
modes  in  Table  2  and  Table  3  for  two-dimensional  and  three-dimensional  cases, 
respectively. 


Table  2  Values  of  Cnm  for  Two-Dimensional  Chamber 


Cnm 

m  =  1 

m  =  2 

m  =  3 

n  —  1 

0.386 

0.0 

0.204 

n  =  2 

0.0 

0.363 

0.0 

n  =  3 

0.204 

0.0 

0.286 

Table  3  Values  of  Cnm  for  Three-Dimensional  Chamber 


Cnm 

m  =  1 

m  =  2 

m  =  3 

n  =  1 

0.72 

0.78 

-0.13 

n  =  2 

0.78 

1.43 

0.0 

n  =  3 

-0.13 

0.0 

0.37 

Using  these  coefficients,  the  direct  numerical  integration  of  Eq.  (5.13)  for  the  first 
two  modes  is  performed  for  the  investigation  of  limit  cycle  behavior  in  a  system. 

A  limit  cycle  behavior  is  found  to  exist  even  though  there  is  complicated 
coupling  among  modes  in  the  baffled  system.  The  result  for  two-dimensional 
rectangular  chamber  is  shown  in  Fig.  5.2,  for  the  first  two  transverse  modes 
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combination;  note  that  the  radian  frequency  of  each  mode  is  normalized  with 
respect  to  the  fundamental  mode  of  an  unbaffled  chamber.  For  a  system  with 
a  positive  growth  constant  of  the  first  mode,  the  stability  of  the  system  can  be 
improved  if  the  growth  constant  of  the  second  mode  is  negative.  Although  this 
condition  corresponds  to  a  decrease  of  the  limit  cycle  amplitude  of  the  second  mode, 
it  may  result  in  an  amplitude  increase  of  the  first  mode,  since  a  higher  first  mode 
amplitude  would  be  required  to  balance  the  stronger  damping  effect  of  the  second 
mode  and  thus  produce  the  stable  limit  cycle.  It  is  interesting  to  mention  that 
the  amplitude  histories  of  limit  cycles  observed  in  this  study  have  different  shapes 
from  those  for  unbaffled  chambers.  The  limit  cycle  amplitudes  are  fluctuating  with 
respect  to  constant  values,  possibly  due  to  the  strong  coupling  between  modes. 

From  energy  considerations,  a  mode  with  a  positive  growth  constant  can  be 
directly  interpreted  as  unstable  if  there  is  no  nonlinear  coupling.  Therefore,  such  a 
mode  can  serve  as  an  effective  source  of  energy  for  the  nonlinear  acoustic  system. 
In  a  similar  manner,  a  mode  with  a  negative  growth  constant  can  be  regarded  as  a 
sink  of  energy.  The  nonlinear  coupling  allows  for  inter-modal  energy  exchange  and 
can  thus  channel  net  energy  from  the  unstable  to  the  stable  modes.  When  the  limit 
cycle  is  reached,  the  sink  and  source  effects  of  the  individual  mode  balance  each 
other  and  the  energy  of  the  system  remains  constant  in  time.  For  the  two-mode 
combination  considered  here,  energy  is  supplied  through  the  first  mode  by  linear 
processes  because  its  growth  constant  is  positive.  This  energy  is  then  transferred  to 
the  second  mode  through  the  second  order  nonlinear  gasdynamics  coupling,  where 
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Figure  5.2  An  Example  Showing  the  Existence  of  Limit  Cycle 
for  Two-Dimensional  Chamber; 
the  First  and  Second  Transverse  Modes 
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it  is  dissipated  by  a  linear  decay  process  acting  on  the  second  mode  with  its  negative 
growth  constant.  In  general,  nonlinear  processes  associated  with  gasdynamics  have 
a  natural  tendency  to  cause  energy  flow  from  lower  to  higher  modes. 

For  a  three-dimensional  cylindrical  baffled  combustion  chamber,  only  three 
cases  of  nonlinear  mode  combinations  are  treated  here.  Investigated  are:  (1) 
the  first  tangential  and  the  first  radial  modes,  (2)  the  first  tangential  and  the 
second  tangential  modes  and  (3)  the  first  tangential,  the  first  radial  and  the  second 
tangential  modes.  These  three  situations  represent  the  most  commonly  observed 
modes  in  actual  rocket  engines,  without  presenting  a  formidable  calculational  task. 

Limit  cycles  are  reached  for  all  cases  and  the  results  are  shown  in  Fig.  5. 3-5. 5 
for  1T/1R  mode,  1T/2T  mode  and  1T/1R/2T  mode  combinations,  respectively. 
For  all  three  cases,  the  limiting  amplitudes  fluctuate  in  a  manner  similar  to  the 
results  for  the  two-dimensional  baffled  chamber,  due  to  the  complicated  coupling 
among  modes.  This  fluctuation  can  be  equivalently  stated  as  a  vertical  shift  of  the 
envelope  of  the  time-varying  amplitude  as  seen  in  Fig.  5.6,  thus  making  the  limit 
cycle  fluctuate  with  respect  to  constant  value. 

The  overall  trends  axe  the  same  as  the  two-dimensional  case,  but  it  is  interesting 
to  observe  that  the  intensity  of  the  fluctuation  for  the  three-dimensional  case  is 
stronger.  The  fluctuation  is  more  than  30  percent  of  the  average  amplitude  for 
all  three-dimensional  cases.  This  means  a  larger  shifting  of  the  envelopes  of  the 
time-varying  amplitudes.  Coupling  among  modes  allows  a  cascade  of  energy  from 
lower  modes  to  higher  modes,  as  evidenced  by  the  limit  cycle  growth  for  the  first 
tangential  mode.  Figure  5.5  shows  that  the  amplitude  of  the  first  radial  mode  is 


PRESSURE  AMPLITUDE 


PRESSURE  AMPLITUDE 


Figure  5.6  Time  History  of  the  Pressure  Amplitude 
for  Three-Dimensional  Chamber;  the  First 
Tangential/Second  Tangential/First  Radial  Modes 
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similar  to  that  of  the  second  tangential  mode.  Thus,  both  the  first  tangential  and 
the  first  radial  modes  are  important  in  dissipating  the  energy  of  the  first  tangential 
mode  for  the  creation  of  a  stable  limit  cycle  in  this  scenario. 

In  summary,  the  typical  nonlinear  behavic.s-limit  cycles  are  observed  in 
combustion  chambers  even  though  there  are  strong  couplings  among  modes  in 
baffled  environments.  And  energy  cascade  through  the  nonlinear  gasdynamics 
processes  is  also  predicted  for  several  mode  combinations.  But  unlike  the  unbaffled 
case,  the  limit  cycle  amplitudes  show  the  fluctuating  characteristics  due  to  the 
complicated  inter-modal  coupling  in  the  baffled  system. 


CHAPTER  6 


CONCLUSIONS  AND  FUTURE  WORK 


Linear  and  nonlinear  acoustic  wave  analyses  have  been  developed  for  two-  and 
three-dimensional  baffled  combustion  chambers.  A  wave  equation  governing  flow 
oscillations  is  derived  from  the  conservation  equations  for  a  two-phase  mixture. 
To  solve  this  wave  equation  in  baffled  environments,  the  eigenfunction  expansion 
technique  is  used  to  express  solutions  in  the  baffle  compartments  and  the  main 
chamber,  separately.  With  appropriate  matching  of  the  acoustic  fields  at  the 
interface  between  the  baffle  compartments  and  the  main  chamber,  unsteady  flow 
structures  in  the  entire  chamber  can  be  obtained.  The  major  advantage  of  this 
approach  is  that  it  alleviates  the  computational  burden  associated  with  conventional 
numerical  methods,  and  provides  more  sweeping  results  in  a  much  more  efficient 
manner.  Furthermore,  it  provides  an  analytical  framework  for  studying  mechanisms 
proposed  as  the  causes  of  unsteady  motions  in  baffled  rocket  engines. 

From  the  linear  acoustic  analysis,  the  detailed  wave  structures  inside  the  baffled 
combustion  chambers,  including  distributions  of  acoustic  pressure  and  velocity,  have 
been  investigated.  Three  stabilizing  effects  of  the  baffle  on  transverse  modes  of 
instability  have  been  observed.  First,  transverse  waves  can  be  longitudinalized 
inside  baffle  compartments.  This  effect  is  more  conspicuous  for  chambers  with 
longer  baffle  lengths.  Longitudinal  waves  are  rarely  a  problem  in  most  liquid 
rocket  engines,  so  the  longitudinalization  of  the  waves  seems  to  be  a  sensible 
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way  for  causing  stability.  Physically,  this  longitudinalization  may  act  to  decouple 
combustion  from  oscillatory  effects  if  the  processes  near  the  injector  face  are 
sensitive  to  transverse  variations  in  pressure.  However,  there  is  also  a  potentially 
destabilizing  concentration  of  pressure  amplitude  near  the  injector  face  associated 
with  longitudinalization  of  the  wave.  Second,  a  severe  restriction  of  the  transverse 
component  of  the  acoustic  velocity  fluctuation  is  observed.  If  processes  sensitive 
to  acoustic  velocity  are  responsible  for  instability,  this  restriction  could  be  the  key 
in  stabilizing  the  system.  Third,  the  normal  mode  frequency  is  decreased  with  the 
addition  of  baffles.  Acoustic  pressure  oscillations  that  might  have  induced  instability 
without  the  baffles  might  fail  to  interact  constructively  due  to  the  shift  of  frequency. 

From  the  nonlinear  spinning  mode  analysis  in  a  cylindrical  unbaffled  chamber, 
linear  couplings  between  modes  are  found  to  exist  and  to  play  an  important  role  on 
the  spinning  transverse  mode  of  oscillations.  Limit  cycle  histories  are  also  calculated 
and  presented.  Nonlinear  acoustic  solutions  for  baffled  combustion  chambers 
are  found  by  using  a  series  of  linear  modes,  with  time- dependent  amplitudes  as 
coefficients  to  express  the  nonlinear  behavior.  A  system  of  second-order  ordinary 
differential  equations  in  terms  of  the  time-varying  amplitudes  is  solved  to  yield 
information  regarding  nonlinear  phenomenon  such  as  limit  cycles,  whose  behaviors 
are  observed  for  several  mode  combinations.  Unlike  the  unbaffled  chamber,  however, 
the  envelope  of  the  time-varying  amplitude  is  shifted,  causing  fluctuations  of  the 
limi'  cycle  due  to  the  strong  inter-modal  coupling  in  the  baffled  system. 


While  this  research  offers  investigation  of  the  detailed  wave  structures,  and 
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some  nonlinear  behaviors  within  the  baffled  combustion  chambers,  some  aspects 
require  further  attention  in  future  research.  First  of  all,  the  mean  flow  Mach  number 
is  considered  to  be  constant  throughout  the  chamber  since  concentrated  combustion 
is  assumed  in  this  analysis.  Consideration  of  distributed  combustion  may  provide 
a  more  accurate  and  realistic  solution  for  this  problem.  Second,  the  study  of 
vortical  structures  by  Wicker32  shows  that  considerable  amounts  of  vorticity  can  be 
generated  by  the  acoustic  wave  from  the  baffle  blade.  Viscous  energy  dissipation  and 
kinetic  energy  convection  by  these  vortices  are  plausible  explanations  of  the  other 
stabilizing  effects  of  baffles  and  are  worthy  of  continued  pursuit.  Third,  flow  turning 
effects  may  be  important  in  the  baffled  chamber  analysis,  since  this  is  a  mechanism 
of  acoustic  energy  loss  due  to  the  turning  of  flow  by  the  acoustic  wave.  Placing 
baffle  blades  inside  the  chamber  can  change  the  direction  of  the  flow  to  generate 
energy  loss  in  a  chamber.  Fourth,  present  nonlinear  analysis  considers  only  the 
second-order  nonlinear  gasdynamics  terms  in  the  formulation.  An  extension  to  the 
higher-order  terms  may  be  able  to  handle  more  realistic  situations.  The  further 
consideration  of  these  four  points  can  aid  in  the  construction  of  a  more  realistic  and 
complete  model  of  actual  baffled  systems. 
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APPENDIX 


CHARACTERISTIC  EQUATION  FOR  THE  FREQUENCY 

InEq.  (3.48),  coefficients  qmiA ,  Cmb,  ?nc>1 ,  qnc  2  and  Cnc  are  all  functions 

of  frequency.  Therefore,  Eq.  (3.48)  can  berewritten  as 


Bn  =  E  G»m-(n)Bm> 


(Al) 


where 


c^{^y)dy 

N 

‘•-os(!Efry)  cos(*ffy)dy' 

X 

If  M  +  1  terms  are  expanded  in  Eq.  (A.l), 

n  =  0  Bo  =  Goo  Bo  +  GoiBi  + . +  GomBm 

n  —  1  B\  =  Gio-Bj  +  Gii  B\  + . +  G\mBm 


(A.  2) 


n  =  M  Bm  =  GmoBm  +  Gm\Bm  + 
In  matrix  form, 


+  GmmBm 


Goo  —  1 

Go\ 

G  02 

Gom 

So 

Guo 

Gu  —  1 

G12 

Gim 

s, 

G20 

G21 

G22  —  1 

Gom 

S2 

Gmo 

Gmi 

G\12 

•  •  •  Ga/a*  —  1  . 

Sa/ 
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That  is, 

[G]  [B]  =  0  (A3) 

For  the  non-trivial  solution,  the  determinant  of  the  matrix  G  should  be  zero. 


G|  =  0 


(A4) 


All  elements  of  matrix  G  are  functions  of  frequency.  Therefore,  frequency  can  be 
determined  from  the  frequency  characteristic  equation  Eq.  (A. 4). 


